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A systematic analysis of matched layers is undertaken with special attention to better understand the remarkable 
method of Berenger. We prove that the Berenger and closely related layers define well posed transmission prob- 
lems in great generality. When the Berenger method or one of its close relatives is well posed, perfect matching 
is proved. The proofs use the energy method, Fourier-Laplace transform, and real coordinate changes for Laplace 
transformed equations. It is proved that the loss of derivatives associated with the Berenger method does not 
occur for elliptic generators. More generally, an essentially necessary and sufficient condition for loss of derivatives 
in Bcrenger's method is proved. The sufficiency relies on the energy method with pseudodiffercntial multiplier. 
Amplifying and nonamplifying layers are identified by a geometric optics computation. Among the various flavors 
of Berenger's algorithm for Maxwell's equations our favorite choice leads to a strongly well posed augmented 
system and is both perfect and nonamplifying in great generality. We construct by an extrapolation argument an 
alternative matched layer method which preserves the strong hyperbolicity of the original problem and though 
not perfectly matched has leading reflection coefficient equal to zero at all angles of incidence. Open problems are 
indicated throughout. 
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1. Introduction 

This paper analyses absorbing layer methods for calculating approximations to the solution, U, of first 
order systems of hyperbolic partial differential equations, 

d 

L{dt,d,)U := dtU + Y^AidiU ^ F, {t,x) G IRl+^ C/(t,a;)eC^. (1.1) 

Approximate values are sought on a finite domain. The source term F and/or initial condition is compactly 
supported in the domain. The absorbing layer strategy surrounds the domain with a layer of finite 
thickness intended to be absorbing and weakly reflective. 

The simplest case is dimension d = 1 with computational domain xi < and absorbing layer in 
xi > 0. For the flrst example consider inhomogeneous initial data and zero right hand side. The simplest 
absorbing layers add a lower order term alxj^>o C U where 1 denotes the characteristic function, for 
example, 

To get a feeling for the reflections, consider the solution U{t,xi) so that, 

for i < , U = {5{xi - t) , 0) . 

Then 

C/i =(5(xi -Oe-""^^-i>«, {dt - di + bal^,>o )U2 = -a l,,>o c?7i. 
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If c ^ 0, then ^t,xiU2 is discontinuous across the ray {xi = —t}. From the perspective of a numerical 
method, such a reflected singularity is undesirable. 

The reflected singularity from a discontinuous lower order term is weaker than the singularity of the 
incident wave. For the equation 

dtU + AidiU + al^.^oCU - 0, 

if C is diagonal in a basis diagonalising Ai, the reflections are avoided. The ease of eliminating reflections 
for this problem with d = 1 is deceptive. No such simple remedy exists in dimensions d > 1. For symmetric 
hyperbolic systems Ai ~ Al, it is wise to choose C = C* > so that the absorption term is dissipative 
in the i^(]R'') norm. 

Consider next the wave equation with friction dttv — diiv + 2a la;i>o dtv = written in characteristic 
coordinates (Ui, U2) = {dtV — div, dtv + div) with absorption B = aC: 

dtU + {^Q^^diU + <7lx,>oCU = 0, ^ = (II 

The absorption matrix C is symmetric and nonnegative but does not commute with Ai. It produces 
unacceptably strong reflections. The absorption from Israeli and Orszag [16], dttv—diiv+a{dtv+div) = 0, 
absorbs only rightward waves and corresponds to 

introducing the notation 7r+(j4i) for the spectral projector on the eigenspace corresponding to strictly 
positive eigenvalues of Ai. The general nonnegative symmetric choice commuting with Ai is a positive 
multiple of 

C = 7r+(Ai) + j/7r„(Ai), iy>0. (1.2) 

We call these smart layers. They dissipate the norm. As observed by Israeli and Orszag, the numerical 
performance of the smart layers is not as good as one would hope. One reduces reflections by choosing 
<7{x) > vanishing to order fc > at the origin. That reduces the rate of absorption and thereby increases 
the width of the layer required. The leading reflection by such smart layers of incoming wave packets 
of amplitude 0(1) and wavelength e is 0(e'^+^). The leading reflection is linear in a. In section [6l we 
introduce the method of Harmoniously Matched Layers which remove the leading order reflections (at all 
angles of incidence) by an extrapolation. 

Open problem. Repeated extrapolation further reduces the order of reflection. It is easy to program 
and it is possible that an optimization could pay dividends. 

Elaborate absorbing layer strategies, like Berenger's PML introduce operators related to but often 
more complicated than the original operator L. The operators in the absorbing layer and in the domain 
of interest may not be the same. For the case of a layer in {xi > 0}, absorbing layer algorithms solve a 
transmission problem for an unknown (V, W) where ^ is a valued function on xi < and is a 
function on xi > 0. The equations in a;i > are chosen to be absorbing and the transmission problem 
weakly reflective. The ingenious innovation of Bcrenger was to realize that the operator R in the layer 
can differ substantially from L. He increased the number of unknown functions in the layer. So W is C^^ 
valued with M > N. 

The pair (V, W) is determined by a well posed transmission problem, 

LV = F on := {{t,x) : xi < 0} RW = on M^+'' , (1.3) 

with the homogeneous transmission condition 



{V,W) e Af on {xi = 0}. (1.4) 
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Here Af C x C*^ is a linear subspaceB The choice of the hyperbohc operator R and transmission 
condition Af is made with three goals, 

• The transmission problem is well posed, and not hard to approximate numerically. 

• Waves from the left are at most weakly reflected at xi = 0. 

• Waves moving rightward decay rapidly in xi > so that the layer can be chosen thin. 

The criterion for perfection that we adopt is that of Appclo, Hagstrom and Kreiss [3] . In the case of 
one absorption, it is formulated as follows. 

Definition 1.1. A well posed transmission problem is perfectly matched when for all F supported in 
xi < , t > 0, the solution supported in t > satisfies V = U\^^^^. 

We prove in ij4. 1.41 that Bcrcnger's method with one discontinuous absorption ai is perfect in this sense. 

In practice one does not absorb in only one direction and the computational domain is rectangular. 
We give in a definition with absorptions in more than one direction and a proof of perfection. 

The strategy of Bcrcnger is quite ingenious. For an artificial boundary in two dimensions at {xi = r} 
and domain of interest {xi < r} it consists of two steps. The first is a doubling of the system and 
the second is insertion of an absorption term in {xi > r}. The doubled system involves the unknown 
U := (C/\ C/^) e C'^ X C^. When F = 0, the doubled equation without dissipation is 

dtW + Aj9j(C/i + C/2) = 0, j = l,2. 

The system with damping in xi changes the j ~ I equation to 

dtU^ + AidiiU^ + U^) + a{xi)U^ = 0, supper C {xi > r} . 

Then U := J^j U'' satisfies L{d)U = in < r. In practice it is the restriction of [/ to xi < r that is 
of interest. There are three distinct ways to view this. One can think of the unknowns as U defined in 
xi < r and U in xi > r with the transmission condition that AiU ~ Ai{U^ + U'^) on xi = r. One is 
given initial values of U and takes initial values of U vanishing. This is the most natural choice and the 
one presented by Bcrcnger. 

From the computational point of view it is simpler to have the same unknowns throughout. The 
simplification is greater when one passes from the half space case to a computational domain equal to a 
rectangular domain in R'*. One introduces U everywhere with transmission condition [Ai{U^ + L/^)] = 
where [*] denotes the jump at xi = r. The transmission condition is then equivalent to the validity of 
the differential equation satisfied by U in all of M.'^. When one uses U everywhere, the initial values of U 
are taken equal to zero outside the computational domain. The initial values are constrained to satisfy 
U = within the computational domain. The choice is otherwise arbitrary. For the case of the 

doubling above the choice C/^ (0,x) = J7(0,a;)/2 for j = 1,2 is common. 

If the domain of interest is |xi| < r one would choose cr > on jxij > r and vanishing for \xi\ < r. 
The transmission condition is [j4i(C/^ + C/^)] = with the jump at xi = r and also at xi = — r. 

In a rectangular geometry in R'' introduce U := (U^, . . . , If^), where G for 1 < Z < d. Then U 
with values in C^'' is required to satisfy (in the case F = 0), 

d 

{L{dt,d^)U)i := dtU' + AidiC^W) + ai{xi)U' = 0, 1 < I < d. (1.5) 

Each absorption coefficient ai{xi) > depends on only one variable. It is strictly positive between the 
inside rectangle and a larger outside rectangle. In the layer between the rectangles the solution is expected 
to decay. If U solves (|1.5|) . then U = solves (|1.1|) on the set {x : ai{xi) = for 1 < ^ < d} 

including the inner rectangle. In the case considered by Bcrcnger the a were discontinuous and the 



''Transmission conditions which involve derivatives can also be treated. The algorithms of Bcrcnger and our HML do not 
require that generality. 
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equations (|1.5p are equivalent to transmission problems where on the discontinuity surface of aj one 
imposes the transmission condition of continuity of Aj U^. 

Our first technique is the energy method. In ^3.21 we show that if (^i , . . . , ^d) = docs not meet 
the characteristic variety of L then the Berenger method is well posed without loss of derivatives. This 
applies in particular to linearized elasticity and suggests that in some ways the Berenger method is 
better adapted to that situation than the Maxwell equations for which it was intended. In §3. 31 we give a 
nontrivial extension of the method of Metral and Vacus to show that Berenger's method for the Maxwell 
equations in dimension d = 2 (resp. d = 3) is well posed provided that (Tj{xj) G W^'°°{Rxj) (resp. 
cTjixj) G W'^'°°{M.x^)). The method introduces a norm that is the sum of L'^{]&f^) norms of suitable 
differential operators Pa{D) applied to U. It has the property that the norm at time ti is estimated 
in terms of the norm at time t2. If one introduces the vector of unknowns Pa{D)U this shows that 
the Berenger problem becomes strongly well posed without loss of derivatives. Such transformations are 
typical of weakly well posed problems, (see the Dominies' proof of Theorem 1.1 in §IV.l of [30]). 

When such an estimate is known, we prove sharp finite speed in §3.41 and perfection in §3.51 and §3.61 
the latter concerned with several variants of the Berenger strategy. The perfection proof passes by a study 
of the Laplace transform on {Imr = 0}. The transformed problem is conjugated to the problem without 
absorption by a r dependent change of independent variable x, an idea inspired by jlOj . 

Our second method is the Fourier-Laplace method. Berenger introduced his PML for Maxwell's equa- 
tions with a piecewise constant. Using a computation which resembles plane wave analysis of reflections 
for problems without lower order terms, Berenger argued that the layers were perfectly matched for all 
wave numbers and all angles of incidence. Using variants of the same approach other closely related PML 
were constructed afterward. Performance is observed to be enhanced using a which are not discontinuous. 
Twice diffcrentiable cubic functions arc the most common. The Berenger method is a very good method 
for Maxwell's equations. The Fourier-Laplace method gives a framework for understanding the compu- 
tations of Berenger. In addition, it is the only method we know for proving well posedncss of Berenger's 
PML with discontinuous a for Maxwell's equations. 

Plane wave analysis is sufficient to study reflection and transmission for linear constant coefficient 
operators without lower order terms. Problems with lower order terms require other tools as it is no 
longer true that the plane waves generate all solutions. The first level of generalization is to use the 
Fourier-Laplace transform for problems where an absorbing layer occupies xi > and both L and R have 
constant coefficients. Hersh |11| found necessary and sufficient conditions for (weak) well posedncss of 
transmission problems. We recall those ideas in ij4.1.1l including the modifications needed for characteristic 
interfaces, and verify in ^4. 1.31 that the condition is satisfied for the Berenger splitting of general systems 
with one discontinuous absorption coefficient. To our knowledge this is the first proof that the Berenger 
split transmission problem with discontinuous a{xi) is well posed. 

We give necessary and sufficient conditions for perfection at a planar interface. In §4.1.41 we verity 
that the condition is satisfied for the Berenger splitting. 

In i 34.1.5l wc prove that in the case of Maxwell's equations (and not in general) the perfection criterion 
follows by analytic continuation from the plane wave identities established by Berenger. 

In ^ 34. 21 we prove using the Fourier-Laplace method that Berenger's method with one coefficient 
cri(a;i) S Lip(R2;j) is well posed and perfectly matched. In our use of the Fourier-Laplace method, 
including this one, a central role is played by the Seidenberg-Tarski Theorem estimating the asymptotic 
behavior of functions defined by real polynomial equations and inequalities. The Fourier-Laplace method 
is limited to coefficients that depend only on xi. 

Our third method of analysis is to study the behavior of short wavelength asymptotic solutions. For 
such solutions we examine in ^JS] the decay in the absorbing layers, and reflections at discontinuities of 
aj{xj) or its derivatives when smoother transitions are used. For problems other than Maxwell, Hu |14| 
and Becache, Fauqueux , and Joly [S] have already shown that the supposedly absorbing layers may in 
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fact lead to growth. The study of short wavelength solutions in the layer yields precise and clear criteria, 
also valid for variable coefficients, explaining the phenomenon. 

The analysis of the reflection of short wavelength wave packets at the interface with the layer also 
leads us to propose in ^J6l a new absorbing layer strategy which we call Harmoniously Matched Layers. 
The method starts with a smart layer for a symmetric hyperbolic system. Then for wavelength e asymp- 
totic solutions of amplitude 0(1) and discontinuous cr, the leading order reflected wave at nonnormal 
incidence typically has amplitude proportional to ae. Therefore an extrapolation using computations 
with two values of a eliminates the reflections proportional to cr. This yields a method with leading order 
reflection 0{e^) at all angles of incidence. The resulting method inherits the simple estimates of the 
symmetric systems. More generally if the first discontinuous derivative of the absorption coefficient is the 
fc*^ then the reflection is 0{[D'^a]e''^^) and the same extrapolation removes the leading order reflection. 
In §6.41 we investigate several implementations of this idea and show that the method with cubic a is 
competitive with that of Berenger with the same a. On short wavelengths or random data it performs 
better than the Berenger method. On long wavelengths Berenger performs better. 

Though we provide satisfactory answers to a wide range of questions about absorbing layers, there is 
a notable gap. 

Open problem. For the original strategy of Berenger for Maxwell's equations with discontinuous ab- 
sorptions in more than one direction we do not know if the resulting problem is well posed. Discussion. 
1. In ij3.6l wc prove well posedness and perfection for a closely related method. 2. In practice discontinu- 
ous a have been abandoned, but it is striking that this problem remains open. 3. Once well posedness is 
proved, perfection follows by the proof in H3.5I 



2. Well posed first order Cauchy problems 
2.1. Basic definitions 

Consider a first order system of partial differential equations for valued functions on M^+'', 

d 

£{x,dt,d,,)U := dtU + ^y^/a/C/ + B{x)U = 0. (2.1) 

1=1 

The principal part of £, denoted Ci, 

d 

Ciidud^) := 9t + 

1=1 

has constant matrix coefficients Ai- In the Berenger strategy, the operators L are the centerpieces and 
they differ from L. It is for this reason that we introduce £ that can be L or L. 

Definition 2.1. The characteristic variety Char(£) C C^^'^ \ {0} of C is the set of (r, ^) such that 
dct£i(r,0 = 0. 

Definition 2.2. The smooth variety hypothesis is satisfied at (t, ^) G Char(£) if there is a conic 
neighborhood 51 of ^ e M'* \ {0} and a C°° function ^ i-^ r(^) on f2 so that on a neighbourhood of (r, ^), 
the characteristic variety has equation r ~ t(^). At such a point the associated group velocity is defined 
to be V := -V^t{^). 

Example 2.1. This hypothesis holds if an only if for ^ near ^ the spectrum of C{0,^) near — r consists 
of a single point with multiplicity independent of ^. For the polynomial (r -I- £,i){t'^ — |^P) with d > 1 the 
hypothesis fails at and only at r + =0 where two sheets of the variety are tangent. Replacing the first 
factor by t + c with c > 1 the hypothesis fails where the two sheets cross transversally. For < c < 1 
the hypothesis holds everywhere. 
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The Cauchy problem for C is to find a solution U defined on [0, oo[xM'^ satisfying (|2.ip with prescribed 
initial data U{0,-). 

Definition 2.3. The Cauchy problem for C is weakly well posed if there exist q > 0, K > 0, and 

a G K so that for any initial values in H'^{R'^), there is a unique solution U € C°([0, +oo[; with 

yt>0, ||[/(^,.)||l2(e.) </^e"*||[/(0,-)||if,(E^). (2.2) 
When the conclusion holds with q ~ 0, the Cauchy problem is called strongly well posed. 

Theorem 2.1. 

(i) The Cauchy problem for Ci is weakly well posed if and only if for each ^ € M.'^, the eigenvalues of 
£i(0, ^) are real. 

(ii) The Cauchy problem for Ci is strongly well posed if and only if for each ^ G M"^, the eigenvalues of 
£1(0,^) are real and Ci{0, ^) is uniformly diagonalisable, there is an invertible S{S,) satisfying, 

S'(C)-i/:i(0, e) S{0 - diagonal, S" , 5"^ G L°°(M^) . 

(Hi) If B has constant coefficients, then the Cauchy problem for C is weakly well posed if and only if there 
existes Af > such that for any ^ G , det C{T,i) = =^ |ImT| < M. 

Remark 2.1. 



1. The algebraic conditions in (i) and (Hi) express weak hyperbolicity, in the sense of Carding. The 



necessity of uniform diagonalisability in (ii) expressing strong hyperbolicity is due to Kreiss [8], |19| . 
2. An application of Gronwall's inequality shows that if Ci satisfies the condition of Theorem 12. 1[ (ii), 
then for all B{x) G L°°{M.'^ ; Hom(C^)), the Cauchy problem for £1 + B is strongly well posed. 



3. By property (ii) if £ is strongly hyperbolic, then every eigenvalue — r of £i(0,.J) is semi-simple. Equiv- 



alently, for any (r, ^) G Char(£) the eigenvalue of Ci{t,^) is semi-simple, i.e. its geometric multiphcity 
is equal to its algebraic multiplicity. It is equivalent to saying that Ker£i(T, ^) = Ker (£i(r, ^))^, or that 
= Ker £1 (t, ® Range £1 (t, ^ . 



2.2. Characteristic variety and projectors for Berenger's L 

To study the Cauchy problem for Berenger's split operators L one starts with a study of the characteristic 
variety. The coefficients of Berenger's operator L are the dN x dN matrices, 



A, 



Vo 



0/ 



Bix) 



(cji{xi)In 



\ ... CFd{Xd)lN/ 



The principal symbol of L is 



f^iAi+TlN ^lAi 

6^2 6^2 + t/jv 



V UAd 



^Ad 



iiAi \ 
C2A2 



^dAd + tInJ 



(2.3) 



Theorem 2.2. 

(i) The characteristic polynomial of L is 

detZi(T,0 



r^('^-i) det L(r,0. 



(2.4) 
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The polynomial associated to the full symbol including the absorption is 



det L(t, = det l( [](r + a^) , 6 + '^j) ' ^2 Y[ir + (^j) , ■ ■ ■ , ^dJlir + a, 

i=l J^d 

If (t, ^) G Chari wii/i r 7^ 0, the following properties hold, 
(ii) The mapping 



(2.5) 



is a linear bijection from KerLi(r, ^) onto KerLi(T, ^). Its inverse is given by 



(Hi) The kernel of the adjoint Li(r, ^)* is equal to the set of vectors $ = ($,■■• , $) such that 
$ e KerLi(T,^)*. T/ie range of Li{t,^) is equal to the set of vectors = (^E'l,-- - , ^'d) swc/i 
i/iai = /or * e Ker Li(t,0*- 

('iuj If moreover the eigenvalue o/Li(r,^) wzi/i r 7^ is semi-simple, the eigenvalue o/Li(r,^) is 
semi-simple. Equivalently, 

KerZi(r,e) e Range Zi(t,0 = C''^ . 
Proof, (i) Adding the sum of the other rows to the first row in the determinant of Li(t^S,) yields, 



detLi(T,e) 



L(r,C) 

6^2 6^2 + tIn 
£.dAd 



N 



Substracting the first column from the others yields, 

L(r,0 
6^2 t/w 

detLi(r,e) 



... tIn 



The first result follows. For the second write, 



detL(r,^) 



6^1 + (T + CTi)/Ar 6^1 

6^2 6^2 + (r + CT2)/JV 



For each z divide the i row by t + ai to find. 



detL(r,0 = n + detLi(l,^i-,. 



. . ^d^d + (r + crd)^Ar 



By formula (j2.4p this implies 



detL(T,0 = n(^ + '^j)'^'i^t^(l'^V'- 



which is equivalent to p.Sp . 
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[ii) Suppose that 7^ $ = (<I>i, • • • , $d) € Ker Li{t,C)- Then, for any I, 

d 

T<^i + iiAiY,<^j = 0, (2.6) 

Add to find 

d 

L,{r,0 Y.'^j = 0. 

Therefore the map <i> n- — J^j maps Ker Li(t, ^) to Ker Li(r, ^). 

If — Oi equation (|2.6p imphes that ah the $j vanish since t 7^ 0. Therefore the mapping is 

injective. 

Let $ G KerLi(0,O- Define 

$j = (2.7) 

This defines an element <I> = ($1, • • • , ^d) in Kerii(0, (,) with 5$ = $, so the mapping is surjective with 
inverse given by (|2.7p . 

(im) Since Li(r, ^)*<f> = (Li(r, ^)*$, . . . ,ii(T, ^)*$) it foUows that the set of $ is inchided in the kernel. 
Since the matrices are square, KerLi(T,^) and KerLi(T,^)* have the same dimension. The set of $ has 
dimension equal to this common dimension which by (ii) is equal to the dimension of Ker Li(r, ^) proving 
that they exhaust the kernel. The last property follows directly from the fact that Range Li(r,^) is the 
orthogonal of Ker Li (t, ^) * . 

(iv) It suffices to show that the intersection of these spaces consists of the zero vector. Equivalently, it 
suffices to show that there is no $ 7^ in KerLi(r, ^) such that 

V* e KerLi(T,0*, (^^Aj$,*) = 0. 

The quantity above is equal to — (^j^*), and $ would belong to (Kcrii(T, ^)*)-'- = Rangeii(T, ^). Since 
T 7^ and Ker_Li(r, ^) n Range Li(r,^) ~ 0, this would imply that $ = 0, leading to a contradiction. □ 

Denote by I1l{t,£_) (resp. IIj^{t,£_)) the spectral projector onto the kernel of Li(t, ^) (resp. Li(t, ^)) 
along its range. For L it is given by 



nL(r,0 = T^/ {zI-L,{t,0) ' dz 



27r« J\z\=p 

with p small. Like the characteristic variety, 11^ depends only on the principal symbol Li. It is charac- 
terized by, 

ni=nL, UlLi{t,^)^0, Li(T,OnL=0, rankn^ = dimKerLi(T,0, (2.8) 

where the r, ^ dependence of 11^ is suppressed for ease of reading. The first three conditions assert that 
nL(T, ^) is a projector annihilating RangeLi(T, ^) and projecting onto a subspace of Kerii(T, ^). That 
it maps onto the kernel is implied by the last equality. 

Proposition 2.1. The matrix IIj^{t, is given by 

nz(r,0 



/'«'^'n.(.,f).. 

T 


T 


TIl{t 






^dAd 

T 


IIl{t 
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Proof. Call the matrix on the right M(t, ^). The properties of the projectors associated to L yield 
formulas for the block of the products 

{M{T,OLi{T,Okj = -^UlL, = 0, {L,{T,OM{T,Okj=-^L,nLLi = Q, and, 

This proves the first three equalities of (|2.8p . Since M projects onto a subspace of KerLi, rankM < 
dimKer Li. Apply M to a vector (^f, 0, . . . , 0) and compare with part (ii) of Theorem 12. 21 to see that the 
range of M contains KerLi(r, ^) so rankM > dimKerLi. This proves the last equality of (|2.8p . □ 

Remark 2.2. 

1. The characteristic varieties of L and L arc identical in r ^ 0. 

2. In particular, the smooth variety hypothesis is satisfied at (r, ^) with r 7^ for one if and only if it 
holds for both, and the varieties have the same equations and the same group velocities. 

3. When the smooth variety hypothesis is satisfied, the spectral projection IIj^{t{£^), ^) is analytic in ^, 
hence of constant rank. It follows that is a semi-simple eigenvalue of _L(t(^), ^) on a conic neighborhood 

ofe 

If the eigenvalue of Li(r, ^) is semi-simple, the kernel and the range of Li{t,^) are complemen- 
tary subspaces as mentioned in Remark [2.11 3.. and the partial inverse Ql{t,£,) of Li(r, ^) is uniquely 
determined by 

QL(T,OnL(r,0 = 0, QL(r,Oii(T,e) = I-ULir^O- (2.9) 
The partial inverse Qi(t,£^) is defined in the same way from Li(r, ^). 



2.3. The Cauchy problem for Berenger's split operators 

Part (z) of Theorem 12.21 proves the following. 

Corollary 2.1. // the Cauchy problem for L is weakly well posed, then so is the Cauchy problem for the 
principal part Li. 

An important observation is that though the Cauchy problem for Li is at least weakly well posed, the 
root r = is for all ^ a multiple root. When there are such multiple roots it is possible that order zero 
perturbations of Li may lead to ill posed Cauchy problems. The next example shows that this phenomenon 
occurs for the Berenger split operators with constant absorption aj. Theorem 12 .41 shows that when r = 
is a root of constant multiplicity of det Li (r, ^) = 0, the constant coefficient Berenger operators have well 
posed Cauchy problems. Cases where the problem are strongly well posed are identified. In the latter 
cases, arbitrary bounded zero order perturbations do not destroy the strong well posedness. 

Example 2.2. 

1. For L dt + di + 82, det L{t, ^) = r + ^1 + ^2- Therefore t = is a root if and only if ^1+^2= 0. 
The doubled system with absorption <t = 1 in xi is 

detL(T,e,-0 = = (r + e+l)(r-e) + e' = + r - ^ 



The roots of detL(r, ^, — ^) arc t = (-1 ± vTT4^)/2. Taking ^ — > —00 shows that the Cauchy 



problem for L is not weakly well posed by Theorem 12. 1[ (Hi) 
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2. More generally if t + ^i +^2 is a factor of dctLi(T, ^) then for a ^ the operator L is not even weakly 
hyperbolic. In this case (|2.5p implies that + r — ^ is a factor of dot L{t, ^, — ^). 

3. This generalizes to linear hyperbolic factors in arbitrary dimension. 

A key tool is the following special case of Theorem A. 2. 5 in [T5] . 

Theorem 2.3 (Seidenberg-Tarski Theorem). IfQ{p,(), RipX)} and S{p,C) are polynomials with 
real coefficients in the n + 1 real variables [p,C,i, . . . , C,n) and the set 

M{p) {C : i?(p,C)=0, 5(p,C)<0} 

is nonempty when p is sufficiently large, define 

p{p) := sup Q{pX)- 

CeM(p) 

Then either p{p) — +00 for p large, or there are a G Q and A ^ Q so that 

p{p) = Ap'^(l + o(l)), p ^ (X. 

Theorem 2.4. Suppose that t = is an isolated root of constant multiplicity m of detLi(r, ^) = 0. 

(i) If the Cauchy problem for Li is strongly well posed, then for arbitrary constant absorptions aj G C, 

the Cauchy problem for Li + B is weakly well posed, 
(ii) If the Cauchy problem for Li is strongly well posed, and if there is a ^ Q such that KerL(0,^) 7^ 
n Ker^j, then Li(0,^) is not diagonalizable. Therefore the Cauchy problem for L is not strongly 

well posed. 

(Hi) If the Cauchy problem for L is strongly well posed and for all ^, KcrLi(0,^) = H Kcr^j, then 

the Cauchy problem for L is strongly well posed. This condition holds if Li(0, dx) is elliptic, that is 
detLi(0,0 7^ for all real 

Remark 2.3. 

1. Part (i) is a generalisation of results in [15] and Theorem 1 in [5]. In the latter paper. Bccachc et al. 
treated the case N = 2 assuming that the nonzero eigenvalues of Li(0,^) are of multiplicity one. They 
conjectured that the result was true more generally Like them we treat the roots near zero differently 
from those that are far from zero. The treatment of each of these cases is different from theirs. The 
tricky part is the roots near zero. We replace their use of Puiseux series by the related Seidenberg-Tarski 
Theorem O 

2. Arbarbancl and Gottlieb [l] proved (ii) in the special case of Maxwell's equations. The general argument 
below is simpler and yields a necessary and sufficient condition for loss of derivatives when the eigenvalue 
of L(0,^) is of constant multiplicity. 

3. Part (Hi) is new, extending a result in the thesis of S. Petit-Bergez [25] . 
Proof. 

{i) For ^ e M'* \ 0, define for p G M+, 

E{p) := max|lm(T) : detL(T,O=0, ^ G M^, = p^} ■ 

Apply the Seidenberg-Tarski Theorem 12.31 with real variables p, C = (Rc t, Im r, ^) and polynomials 
R{pX) = |detZ(T,^)p + (l^p - p'^y. S ^0 and Q(p,C) = Imr. Conclude that there is an a ^ and a 
rational r so that 

E{p) = a p"" {1 + o{l)) , p^oo. 

To prove the result it suffices to prove that Imr is bounded, i.e. to show that r < 0. Suppose on the 
contrary that r > 0. 
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Given t,^ define k G S'^-\ peR+, and, 9 by 

i r 

Choose sequences T(n), and ^(n) so that for ji oo, 

detZ(r(n),^(n)) = 0, ImT(n) = a(p(n))''(l + o(l)) . (2.10) 

Write 

Z(t,0 = Zi(r,C) + S = p{Liie,k) + ^B^ = p{ei2Nx2N + Li{0,k) + -^By 

The matrix L(t, ^) is singular if and only if —0 is an eigenvalue of Li{0, k) + p^^B. 

For large p this is a small perturbation of Li (0, fc). Choose /.i > so that for |fc| = 1, the only eigenvalue 
of Zi(0, k) in the disk |6'| < 2/z is 6* = 0. 

Because of the strong well posedness of L, there is a uniformly independent basis of unit eigenvectors 
for the eigenvalues of ii(0, fc) in \6\ > p.. By part {iv) of Theorem 12.21 there is a uniformly independent 
basis of unit eigenvectors for the eigenvalues of Li{0, k) in \9\ > p. 

It follows that there is a Cq so that for p > Cq the eigenvalues of Li{Q,k) + p^^B in \9\ > p differ 
from the corresponding eigenvalues of Li(0, fc) by no more than C'o/p. In particular their imaginary parts 
are no larger than Co/p. Therefore, the corresponding eigenvalues t ~ p9 have bounded imaginary parts. 
Thus for n large, E{p{n)) can be reached only for the eigenvalues —9[n) which are perturbations of the 
eigenvalue of Zi(0, k{n)). 

Perturbation by 0(1/ p) of the uniformly bounded family of dN x dN matrices, ii(0, fc), can move 
the eigenvalues by no more than 0{p^'^). Since the unperturbed eigenvalue is 0, |6'(7i)| < C p{n)~^^'^^ , 
so 

|r(n)| < Cpin)^-^, ImT(n) = a p{ny {I + o{l)) , a^O. 
Therefore r <1- 1/dN < 1 and 

d 

l[{r{n) +a,) = r(n)'^(l + o(l)), ^n) l[{r{n) + a,) = ^n) rW^-^ll + o(l)) . 

Insert in identity (|2.5p to find 

detLi(rH'^(l + o(l)), eHT(n)'^-i(l + o(l))) =0. 
Divide the argument by p{n) T{n)'^~^ and use homogeneity to find 

.(^^(l + o(l)),fc(n)(l + o(l))) = 0. 

The constant multiplicity hypothesis shows that 

detLi(r,e) = r"Fi(r,e), and G Fi{0,£_)^0. (2.11) 

Since for n large (r(n)/p(n))(l + o(l)) 7^ we have 



detLiI 



(^^(1 + 0(1)), fc-(n)(l + o(l))) = 0. 



Passing to a subsequence we may suppose that the bounded sequence fc(n) — > fc. In addition, T(n)/ p{n) 
so passing to the limit yields Fi(0, fc) = contradicting p. lip . This contradiction proves (i). 
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(ii) Part (i) of Theorem 12.21 shows that is an eigenvalue of ii(0, ^) with algebraic multiplicity equal to 
A'^(d — 1) + m. It remains to see that with the assumption, the dimension of Ker ii(0, ^) is strictly smaller 
than N{d— 1) + m. By definition 

d 

KerZi(0,C) = {$ = : ^ G HpKer (^p^^)} . 



Define 



£i 



a 



Then £i C Kerii(0,^) and dimfi = N{d- 1) 
Define 



If 



£2 ■■= KerLi(0,C)(»O 
$ e KerZi(0,O, Ej=i e KerLi(0,O , write 



dim £^2 = TO- 



^$,,0,...,0 



e£-2, We£i. 



Thus, KerLi(0,0 C fi ©£:2- 

Pick V in KerLi(0,,^), but not in n Ker^j . Then 

y = (F,0,--- ,0) e ^2 and y^KerZi(0,O- 

This proves that Kerii(0,^) is a proper subset of £1 £^2, so 

dim(KerZi(0,O) < dimfi + dim £2 = N{d-l) + m. 

Thus the geometric multiplicity of the eigenvalue is strictly less than its algebraic multiplicity. There- 
fore, Li(0,^) is not diagonalizable. This proves (ii). 



{Hi) To prove that the split problem is strongly well posed it suffices to consider the principal part. 
Suppose i(0, ^) is uniformly diagonalisable on a conic neighbourhood of ^ G R^\0. For U = {U^, • • • , W^), 
introduce 



V={V^,--- ^V^) with F^:=^C/^ and V' := for 2<l<d. 



(2.12) 



Then 



Li{dt,d,)U = dtV + Q{d^)V ^0, with, QiO:-- 



/ii(0,O o...o\ 

6^2 ... 



V . . . 0/ 



The eigenvalues of Q{£,) are those of Li(0,^), therefore real. It suffices to diagonalize uniformly Q(^) 
on the conic neighbourhood of ^. By homogeneity it suffices to consider ^ with |^| = 1. By hypothesis 
there exist a real diagonal matrix -D(^) and an invertible matrix S'(^) so that 

Li{0,O = S{OD{OS~\O, and, > 0, G ||5(0I1 + H^^'lOII < 

Seek a diagonalization of Q{£,) on |^| = 1 in the form, 



/ sio ... o\ 

^2A2Q{0S{0 Id... 



\UA,QiOS(0 ... IdJ 



so, (sior 



( {sior' ... o\ 

V-^d^<iQ(0 ... idj 



(2.13) 



. . . 0\ 
... 

\UA4l-QiOLi{0,0)S{0 0...0/ 

and 

S-\OQiOSiO isdiagonal ^ - Q(Oii(0, 0) = 0, 2<j<d. (2.14) 

From the strong well posedness of Li it follows that uniformly in ^ one has KerLi(0,^) © 
Range Li(0,^) = C^. Choose Q equal to the left inverse of ii(0,^) defined in Since KerLi(0,0 = 
CiKei ^jAj, the condition on the right in p.l4p holds so diagonalizes Q{£,). Since S'(^) and 5'(^)~^ 
are bounded on a conic neighborhood, it follows that S{^) and S{£,)^^ are bounded on a neighborhood 
of ^ in 1^1 = 1. A finite cover of the sphere, completes the proof. □ 

Remark 2.4. Denote by S{£,) the function homogeneous of degree zero given by (|2.13|) for |^| = 1 with Q 
constructed in the proof. Then \\ S {D)V {t)\\ 1^2 (jg^d-^-^ with V from (|2.12p is a norm equivalent to ||J7(t)||i2(Rd) 
and is conserved for solutions of Li{d)U = 0. Those solutions yield a unitary group with respect to the 
norm \\S{D)V\\L2(^^dy 

3. Analysis of the Berenger's PML by energy methods 

This section contains results proving that the initial value problems so defined are well posed. We begin 
with the case of Gevrey absorptions, then and finally the case of the Heaviside function. 

In Section 13.11 we prove that when L is only weakly well posed, Gevrey regular ui lead to well posed 
initial value problems in Gevrey classes. Commonly used a are not this smooth. 

The strongest result, from Section [3T2l applies when L{Q,dx) is elliptic. Important cases are the wave 
equation and linearized elasticity. In these cases the operator Li is strongly hyperbolic so remains strongly 
hyperbolic even with general bounded zeroth order perturbations. Thus for bounded (Ji{xi) the initial 
value problem is strongly well posed. 

In Section [3731 we analyse the case of L associated to Maxwell's equations with finitely smooth a. We 
follow the lead of [25] and extend the analysis of [22| to several absorptions cr/ and to higher dimensions. 
Related estimates for the linearized Euler equation have been studied by L. Metivier [2T|. 

The results of this section do not treat the case of Berenger's method for Maxwell's equations with 
discontinuous aj. The case of one absorption is treated in §4. A closely related method is treated by an 
energy method in fj321 

3.1. General operators and Gevrey absorption 

The next result is implied by Bronstein's Theorem [6] [7] [23] [24]. It shows that when ii(0,^) has only 
real eigenvalues and the aj belong to the appropriate Gevrey class, then the Cauchy problem for L is 
solvable for Gevrey data. 

Definition 3.1. For 1 < s < 00, / e S'iW^) belongs to the Gevrey class G"'(M'^) when 

3C,M, VaeN'', < Afa!Cl"l. 

Then C f^„H'^{W^) C C°°(R''). For s > 1 the compactly supported elements of G^ are dense. If 
1/(01 < Ce-I^l" with < a < 1, then u e G^'". 

Theorem 3.1. // the principal part Li is weakly hyperbolic, and aj G G^^^'^^(R'') then for arbitrary 
f e G^/(^+i)(M'^) there is one and only one solution u e C°°{m}+'^) to 

Lu = 0, «(0,-) = /. 

The solution depends continuously on f. 
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3.2. Strong hyperbolicity when L{0,d) is elliptic 

Theorem 3.2. // L is strongly well posed and L{0^ d) is elliptic, then L is strongly well posed for any 
absorption {ai{xi), . . . ,(Td{xd)) in (i°°(R))'*. 

Proof. Krciss' theorem 12.11 asserts that an operator with eonstant coefficient principal part is uniformly 
well posed if and only if the principal part is uniformly diagonalisable on a conic neighborhood of each 
^ 7^ 0. Therefore the Corollary follows from the third part of Theorem 12.41 □ 

Example 3.1. This result implies that the PML model for the elastodynamic system is strongly well 
posed. The system is written in the velocity-stress {v, S) formulation, 

pdfV — div E = , 9tS — Ce{v) = , {diVj + djVi) , 

with positive definite elasticity tensor C and E := Ce. See [S], where the authors showed that such layers 
may be amplifying (see Section [5]). 



3.3. The method of Metral-Vacus extended to the 3d PML Maxwell system 

Metral and Vacus proved in [22] a stability estimate for Berenger's two dimensional PML Maxwell system 
with one absorption cti(xi) G W^''°°{M.) and x = {xi,X2) € K^. There are two crucial elements in their 
method. First following Berenger, they do not split all variables in all directions. This section begins 
by showing that the partially split model is equivalent to the fully split model restricted to functions U 
some of whose components vanish. The L evolution leaves this space invariant and its evolution on that 
subspace determines its behavior everywhere. 

The second element is that on the partially split subspace there is an a priori estimate bounding 
the norm at time t by the same norm at time 0. This looks inconsistent with the fact that the Cauchy 
problem is only weakly well posed. However the norm is not homogeneous. Certain linear combinations 
of components have more derivatives estimated than others. The observation of [35] is that the system 
satisfied by the fields and certain combinations of the fields and their derivatives, yields a large but 
symmctrizable first order system. These estimates have been obtained, and extended in Sabrina Petit's 
thesis [25j in the 2d case with two coefficients, and in the 3d case for an absorption in only one direction. 

In this section, motivated in part by the clarification of the role of symmetrizers in the work of L. 
Metivier [3T] for the 2d variable coefficient Euler equations in geophysics, we construct analogous more 
elaborate functionals which suffice for the general case of three absorptions in three dimensions. They 
require a, G W'^-°°{R). 

Maxwell's equations for dtEi and dtBi contain only partial derivatives with respect to X2,X3 and 
not xi. In such a situation Berenger splits the corresponding equations in directions X2,X3 but not in 
direction xi. To see why this is a special case of the general splitting algorithm (|1.5p reason as follows. If 
the equation for dtUj from L does not contain any terms in dk, that is the j"^ row of vanishes, then 
the equation for the j'*^ component of the unknown U'^ corresponding to the splitting for the fc"^ space 
variable is. 



dtU^ + ak{xk)U^ = 0, = e-"'^("'^)*?7fe(o,a;). (3.1) 

Plugging this into the other equations reduces the number of unknowns by one. The simplest strategy 
is to take initial data U^{0,x) = which yields the operator L restricted to the invariant subspace of 
functions so that =0. Conversely if one knows how to solve that restricted system then the full system 
can be reduced to the restricted system with an extra source term from (|3.ip . 

Summary. To study the fully split system it is sufEcicnt to study the system restricted to {Uj = 0}. 
Performing this reduction for each missing spatial derivative, corresponds to splitting equations only 
along directions containing the corresponding spatial derivatives. 
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An extreme case of this reduction occurs if an equation contains no spatial derivatives, that unknown is 
ehminated entirely. For the Maxwell system which is the subject of this section this does not occur. The 
use of unsplit variables 

• reduces the size of U reducing computational cost, 

• corresponds to Berenger's original algorithm, 

• is important for the method of Mctral-Vacus which takes advantage of the vanishing components 

Consider the 3c? Maxwell equations, 

dtE -V X H = 0, dtH + V X E ^0. 
Defining U = E + iH, they take the symmetric hyperbolic form (|l.ip with hermitian matrices 






Ai = , ^2 = and ^3 = i . (3.2) 



Introduce the splitting ()1.5|) with some components unsplit. Define the subspace T-L with vanishing 
components corresponding to the unsplit components 

n := = iU\U^,U^) e H^iR^ ; C^f} : ^ 0, ^ 0, C/| o|. 

For U = {U^,U^, U^) in define 

U := U^ + U^ + U^, := djU, V'^ := dijU, W := J2k<^k{xk)U'' , 

:= d,W, Z := Ekdk{Wk+akixk)Uk), d,Z, (3.3) 



V 



(u, V\V''^ , ,W,W, Z3^ 



-<54 



The function Z and therefore Z^ are C valued. The other slots in V are C"^ valued. The second derivatives 
y*'^ of U are ordered as y^'^, y^'^, V^'^ , V^'^, V^'^, V^'^. This convention is important when the equations 
for V are written in matrix form. Computing in turn , Z, Z^ requires two derivatives of aj . 
The unknown in ^^5^ is C/ = (C/\ C/^ U^). The appear in the fifth slot of V. Therefore, 

•)|1(L2(K3))" > •)il(L2(R3))0- 

For the Cauchy problem the initial data is Uq = (J7g , Uq,Uq), from which Vq is deduced by the derivations 
above, and 

|Vo||(L2(R3))54 < C||C/o||(/i-2(Jj3))3. 

Theorem 3.3. Ifjjj, for j = 1,2,3, belong to W'^^°°{R), then for any Uq = {U^,U^,U^) in U there is 
a unique solution U in L^{0,T;'H) of the split Cauchy problem (|1.5p with initial value Uq. Furthermore 
there is a C'l > independent of Uq so that for all positive time t, 

l|C/(i,-)ll(L2(R3))o < Cie^i*||C/o||(^,(„3))9- (3.4) 

Proof. The main step is to derive a system of equations satisfied by Y{t, x) together with a symmetrizer 
S{D). These imply an estimate for t >0, 

||V(t)||,,(^3) < C2e^^*||V(0)||^,(^3). (3.5) 

From this estimate it easily follows that the Cauchy problem for the V-equations is uniquely solvable. It 
is true but not immediate that if the initial values of V are computed from those of U then the solution 
V comes from a solution U of the Bcrengcr system. The strategy has three steps, 

• Discretizc the Bcrengcr system in x only. 
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• Derive an estimate analogous to p.Sp for the semidiscrete problem. The estimate is uniform as the 
discretization parameter tends to zero. The proof is a semidiscrete analogue of p.4p . 

• Solve the semidiscrete problem and pass to the limit to prove existence. 

This is done for the case d = 2 in [25] to which we refer for details. Uniqueness of the solutions to the 
V-systcm and therefore U is simpler and classical and is also in |25) . 
Equation (jl.5p yields, 

dtW + ^kV^ + <^jU^ = 0- (3-6) 

k 

Summing on j yields 

dtU + LiO,d)U + W ^0. (3.7) 

Differentiate in direction Xj to find, 

dtV^ + L(0, d)V^ + = 0. (3.8) 

Differentiate once more to get 

dtV''^ + L(0, d)V''^ + diW' = 0. (3.9) 
The quantity diW^ on the left is replaced using the next lemma. 
Lemma 3.1. 

d,W ^ L{0, d)A,W + Ze, - Efc E,u{o'^U + guV''), (3.10) 

where Cj is the j^^ vector of the standard basis, and Eij is the 3x3 matrix all of whose entries vanish 
except the {i,j) element that is equal to 1. 

Proof. 

First evaluate L{Q,d)AjW to find 

L{Q,d)A,W = Y,AkA,dkW. 

k 

The matrices in Maxwell's equations satisfy AjA^ = —Ejk for j ^ k, and A'j = J2k^j -^kk = I — Ejj. 
This yields 

L(0, d)AjW = - ^ EjkdkW + (/ - Ejj)djW = djW - ^ EjkdkW = d^W - div{W) e,. 
k^] k 

Introduce the definition of Z to find 

:= djW = L{0, d)AjW + div{W)ej 

= L(0, d)A,W + Ze, - iJ2^ dk{(TkUk)) e.j. 

Compute 

( Efc dkicfkUk)) ej = J2k EjkdkicTkU) 

^T.kE,k^'kU + Y.kE,kakV\ 

which proves (|3.10p . The proof of the lemma is complete. □ 



Differentiate p.lOp in space to obtain 

d,W^ = L(0, d)A,W' + d^Ze, - Y.k Ejkd, {auU + a'^y'^). 
Inserting into (|3.9|) yields 

dtV'^' + L(0, a)r^^ + L(0, d)AjW' + Z'^ - ^ E,kd,{akU + a'kV'') - 0. 
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(3.11) 



This is equivalent to 

dtV'^ + L{0,d)V''^ + L{0,d)AjW' + 

To close the system it remains to evaluate the time derivatives of W, and djZ. 

dtW = Y^akdtU" = -Y,<Jk{AkU + akU") = -Y^akAuV" -Y^alu''. 

Using the particular form of the equations yields 

^ {Y.k^k)W - diag{a2a3,aia3,aia2)U - diag{ai, (72, (T3)W . 

Therefore 

, (Ticra, cri(T2)?7 — diag{(Ji,a2,(J3)W = . (3-12) 

fe k 

Differentiate in Xi to find 

dtW' + Eki^^i'^k)AkV'' + cTkAkV''') + <Jk)W + 'JkW' 

-di{diag{(72(^3,aia3,(7i(72))U - diag{a2(T3,ai(73,aia2)V'- (3.13) 
—di{diag{ai, a2, a3))W — diag{ai, a2, a3)W^ = 



Next compute 



dtZ = J2^^^t{W^ + a,U,). 



Consider the pair of equations 

dtU^ + id2U3 + CJ2UI = 0, and dtUf - id3U2 + aaUf = 0. 
Add the two equations. Also add a2 times the first to (T3 times the second. This yields two equations, 

dtUi + t{d2U3 - 93C/2) + Wi^O, and dtWi + i{cT2d2U3 - as^sC/s) + (rluf + apl = 0. 
Rewrite the last term as (T|C/f + (T|C/f = ((T2 + cra)!^! — (T20'3J7i, to find 

dtUi + i{d2U3 - 831/2) + 14^1 = 0, and, 9(1^1 + i{(J2d2U3 - 0383112) + {^2 + 0-3)^^1 - cT2a3Ui = 0. 
Multiply the first equation by ai and add the second to obtain 

dt{Wi + criC/i) + i ((ffi + a2)d2U3 - (ffl + CJ3)d3U2) + (Efc '^k)Wi - CT2CT3U1 = . 

The other indices follow by permutation. Differentiate in and add to find 

Y^k^kdtiWk +<JkUk) +iJ2k(^k ((o-fc + a-k+i)dk+iUk+2 - (crk + (^k+2)dk+2Uk+i) 

+ ^^ HEk <^k)W,) - Efe dk{ak+i<Jk+2Uk) = 0. 
The terms with two spatial derivatives cancel. This leaves 
dtZ + iY.k o'k{dk+iUk+2 - dk+2Uk+i) 

^k){Y.k w^) + Efe <Wk Ek <^k+iak+2Vt = 0. 

Since 

Z^Y. + <'k{xk)Uk) = Y^W^ + '''kUk + ^kV^). 

k k 

we can replace (Efe o'fc)(Efc ^k) '^^ the previous equation by (E^ tJk){Z - Eki'^'k^k + <7k Vj:)) so 

dtZ + {Ek ^k)z + ^ Ek - Vkti) 



<Ek <^k){Eki<^'kUk + VI:)) + E, <Wk Ek <yk+i<yk+2V, 
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'k+2 



k+2,j 



fe+1 



Replace Z by EA;(^fc + '^kV^) to end up with 

dtz^ + (E. '^^)^^ + Efe(H^ + ^fe^) + ^'^Z' - v^) + ^ E. - ^^') 



Efe ^'kWl - E/c dj{cJk+i(Jk+2)V^ - E/c '9jK+lCrfc+2)^^' 



0. 



Summarizing, V is solution of a first order system, dt¥ + P {dx)y + B {x)Y = 0, whose principal symbol 
is given by 



Pid) 





L{0,d) 


04,6 


® 03,3 


04,36 


§03,3 


04,36 


503,3 


04 


46 


5 03,3\ 


06,4 


® 03.3 


h ® 


L(0,a) (/6®L(0,a))Af 06,3<5 


503,3 


Oe 


46 


503,3 


03,4 


® 03.3 


03,6 


® 03,3 


03,3 <5 


503,3 


03,36 


503,3 


O3 


46 


503,3 


03,4 


«)03,3 


03,6 


«)03,3 


03,3$ 


503,3 


03,36 


503,3 


O3 


46 


503,3 


y 04,4 


'8 03,3 


04,6 


® 03,3 


04,3$ 


503,3 


04,3(5 


503,3 


O4 


46 


5 03,3^ 



M 



Here the V"'^ are ordered as indicated before the theorem and, 

fAi \ 
Ai 



^2 
\0 A3J 

To symmetrize it suffices to construct a symmctrizcr for the upper left hand block 

/h<E)L{0,d) 04,6® 03,3 04,3^03,3 \ 



Qid) 



06,4 «) 03,3 h ® HO, d) {h ® L{0, d))M 
V 03,4 ® 03,3 03,6 65 03,3 03,3 ® 03,3 / 



We verify that 



h^h 04,6® 03,3 04,3^03,3' 
S := I 06,4 «) 03,3 /6®/3 {Ie®l3)M 
\03,4 03,3 03,6 ® 03^3 I3 ® I3 



h^h 04,6® 03,3 04,3® 03,3 

with S'"^ = I 06,4 «) 03,3 /6«'/3 ~M{h®l3) 

\03,4 ® 03,3 03^6 ® 03,3 h ® h 



is a symmctrizcr for Q{i£,)- Compute 
which is symmetric since L{Q, ■) is. 



^/4®L(0, •) 04,6^03,3 04,3® 03,3' 

06,4 ® 03,3 h ® HO, ■) 06,4 «) 03,3 

.03,4® 03^3 03,6 03,3 03,3 (g) 03,3; 
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Therefore P(^) is symmetrizable by a matrix independent of ^. Hence, the Cauchy problem for (|3.7p . 
(|3.8p . (|3.1ip . p.l2p . (|3.13p . (|3.14p is strongly well posed. The norm of the zero order terms depends on 
the coefficients aj and their derivatives up to order 2. The estimate of the Theorem follows. 

Remark 3.1. We recall a computation from [25], showing that when there are only 2 coefficients ci and 
fT2 only one derivative of aj is needed. This is always the case in dimension d = 2. When 0-3 = 0, split W 

as 



W = E33W 




Then 





(<J2Ul 











Rewrite (13.71) as 



= diag{a2, (Ji, 0)U. 



dtU + L{0,d)U + E33W + diag{a2,cri,0)U ^ 0. (3.15) 
Differentiate with respect to xi and X2 to obtain 

dtV^ + L{0, d)V^ + EssdjW + dj{diag{a2,c7i,0))U + diagia2,cri,0)W = 0. (3.16) 
To find an equation on W, proceed as in the 3d proof to get, 

dtW + J2 '^kAkV^ + (^ nk)W - aia2U = (3.17) 

k k 

Therefore V is solution of a first order system, whose principal symbol is given by 



P{d) 



fL{0,d) 
L(0,9) 





V 



A symmetrizer is given by 



/I 













/ 





lE23 








/ 


iEi3 


Vo 








I J 



with 





E33di 
L(0, d) E33d2 

/ 



s- 



/ 1 
0/0 

0/ 

V 



~ 

-iEu 
h . 



3.4. Sharp Finite Speed for Berenger's PML 

Recall some notions associated with estimates on the domains of influence and determinacy for a hy- 
perbolic operator L (see [17j). The timelike cones are the connected components of (1,0, ...,0) in the 
complement of the characteristic variety. The forward propagation cone is dual to the timelike cone. A 
lipschitzean curve [a, 6] 9 f n- {t,^(t) is an influence curve when 7' belongs to the propagation cone for 
Lebesgue almost all < G [a,b]. 

Theorem 3.4. Suppose that L defines a strongly well posed Cauchy problem and that the multiplicity 
0/ r = as a root of det Li{t,^) = is independent of ^ G M"^ \ 0. The support of the solution of the 
Berenger transmission problem is contained in the union of the propagation curves of L starting in the 
support of the source terms when either of the following conditions is satisfied. 

[i) e M"^ \ 0, KerLi(0,O = Ker^^Aj, and Vj, dj e L°° . 

[ii) Li is Maxwell's equation and\fj, aj G W^'°° . 
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Proof. The characteristic varieties of L and L satisfy Chari = CharL U {t = 0}. When {t = 0} has 
muhiphcity as a root of det Li{t,£) = independent of ^ G M'' \ the timehke cones of L and L coincide. 
Therefore the propagation cones and influence curves coincide too. 

Case (i). Part [ii) of Theorem 12.41 proves that Li defines a strongly weU posed Cauchy problem. It 
follows that the sharp propagation conclusion of the Theorem is valid for Li + B{t^ x) for any bounded 
B{t,x). This follows on remarking that the solution of (L + B)U = with initial data Uq is the limit at 
1/ — )• oo of Picard iterates W^. The first, U^, is defined as the solution of the Cauchy problem without B. 
For v > 1 the iterates are defined by, 

Zit/'^+i + B{t,x)U'' = 0, U^+^iQ,-) ^ Uo. 

Since Li has constant coefficients, sharp finite speed is classical for that operator. An induction proves 
that each iterate is supported in the union of influence curves starting in the support of Uq. 

Case (ii). Reason as above constructing by Picard iteration approximations V converging to the 
solution V from p.3p . Since the equation satisfied by V is strongly well posed the iterates converge. An 
induction shows that they are supported in the set of influence curves starting in the support oi Uq. □ 



3.5. Proof of perfection for Berenger's PML by a change of variables 

This section continues the analysis of Berenger's method when the hypotheses of Theorem l3.4l are satisfied. 
In those cases well posedness is proved by an energy method. In addition suppose that 

Vj, 3Lj>0, crj=0 when \xj\ < Lj . (3.18) 

Denote by i? := Ilj] — Lj,Lj[. 

Definition 3.2. In this setting the method is perfectly matched when for arbitrary F G (jO,oo[xR) 
the unique solutions V and U of 

LV = F, Vl^^ = 0, LiU = F, Ul^^ = 0, (3.19) 
with L as in (|1.5p satisfies 

yW.n = Ul.H- (3-20) 

Theorem 3.5. With the assumptions of Theorem \3.4\ and Uj as above, Berenger's method is perfectly 
matched. 

Proof. Taking the Laplace transform of the V equation in (|3.19p yields a transform holomorphic in 
Rer > To with values in L'^(R'^) satisfying for 1 < j < d, 

V'^ + {T + a.i{xj))-^AjdjV ^ f'^, with V:=J2^^^ ^■=J2^^- (^'^l) 

j j 

Multiply by t and sum on j to obtain 

rV> + ^-^^A,a,F = ri?. (3.22) 

When T is fixed real and positive this equation can be transformed to the corresponding equation 
without the aj by a change of variables. The change of variables depends on r. The resulting equation is 
exactly that determining U :— J^j ■ In. this way we find that V is obtained from U by this change of 
variables. This idea is inspired by Diaz and Joly in fTO\ . 

For real r > define d bilipschitzean homeomorphisms Xj{xj) of R to itself by 

dX,{x,) ^ r + a,{x,) ^ X,(0)=0. 

dXj T 
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Then, 

d dXj d T + crj{xj) d T d 



dxj dxj dXj T dXj ' t + aj {xj ) dxj dXj 

Therefore if U{X) is the solution of 

tU{X) + Y,A,^U = FiX), (3.23) 

then the solution V of (|3.22|) is given by V{x) := lJ{X{x)) since the latter function of x satisfies the 
equation determining V . 

Since X{x) = x for x £ R this proves that the transforms of U and V satisfy for real t > tq 

^V^(t,x) = Y.W{t,x), xGR. (3.24) 

j j 

Since both sides of the identity in p.24p are holomorphie in Re r > tq it follows that the identity extends 
to that domain by analytic continuation. 

Equation (|3.2ip and its analogue for U then imply that for all j, V^\^ ~ ^"'Ifl- Uniqueness of the 
Laplace transform implies V-'\j^ = U^j^ for all t proving perfection. □ 

Remark 3.2. The proof is very general. It shows that once the initial value problem defined by L is well 
posed there is perfect matching. The proof works more generally for at least weakly well posed methods 
for which the Laplace transform can be reduced to (|3.22l) for real r. Our favorite version of the Bcrenger 
algorithm is analysed this way in 



3.6. Perfection for methods related to Berenger's PML 

Consider (|3.22l) with F ~ 0. This equation is the starting point for many authors to construct well posed 
PML. It has been viewed as a complex stretching of coordinates (see [29], [^, [26], [H]). This idea, for r 
real, becomes an honest change of variables as in [10] , that is at the heart of the proof in In the ease 
of Maxwell system, it can be viewed as a system with modified constitutive equations (a lossy medium 
[27], [2]), or recovered as above from the Berenger's system. The system p.22p is not differential because 
of the division by r + aj {xj ) . In order to recover a hyperbolic system, a change of unknowns is performed. 
We adopt the approach in [20] for the Maxwell system. 

Lemma 3.2. With matrices given in (j3.2p . define Sj := {t + aj{xj))/T. There exists a pair of invertible 
matrices M, N , unique up to a multiplication by the same constant, such that 

S;^NAj ^ AjM, j = 1,2,3. (3.25) 



They are given by 



'S2S3 

M = 7|GS'2 0|, iV = 7| S1S3 j, 7eC\0. (3.26) 

S1S2/ 




Proof. Since 



Ajej=0, AjCj+i = -iej+2, Ajej+2 = iej+i, 

it is easy to see by applying p.25p to Cj that M is necessarily diagonal, M = diag(TOi, 7712, 'm-s)- Applying 
(|3.25p to Cj^i and ej+2 shows that for any j, 

Ncj+i = mj+2<5'jej+i, Nej+2 = mj+iSjej+2 ■ 

This implies that N is also diagonal, equal to diag(m2S'3, msS*!, TOiS'2), and 

TOi^a = msS*!, m2Si=miS2, 171382=^2283 
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This leaves no choice but to choose (|3.26[> . 
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In the rest of the analysis take 7 = 1. In p.22[) with = replace V replace by U. Insert p.25|) to 
obtain 



tNU + ^Ajil/ajt/ = 0. (3.27) 

j 

The fact that aj depends only on Xj and the form of the matrices guarantees AjdjM = 0. This yields 

AjMdjU = AjdjiMU). 
Define a new unknown V := MU to find 

NM-'^tV + ^AjdjV ^ 0. (3.28) 

j 

NM^^ = dia.g{Si^S2S3, 82^8381, S^^SiS2)- Next compute a rational fraction expansion of T5f ^525*3 
as 

(r + (T2)(t + 0-3) , / , \ , +0-20-3 - (71 (0-2 + 0-3) 

■ = T + (02 + 03-01) + ; . 

T + oi r + oi 

Introduce a new unknown W by 

TNArW = tV + J:W + J:^W, equivalently Wj = ^-^^Vj ^-Uj, 



with 



S := diag(oi,02,03), 

diag(o2 + 03 - 01,03 + oi - 02,01 + 02 - 03), 
-.^ diag((oi - 02)(oi - 03), (02 - oi)(o2 - 03), (03 - oi)(o3 - 02)). 

This leads to a system in the unknowns V and W 

L{dt,d^)V + S^^V + ^ 0, dtW + i:W - V = 0. (3.29) 

Finally, U is recovered from 

U ^ dtW ^ V - T.W . 

The system of equations for V, W is strongly well posed since L is symmetric hyperbolic. In the case of 
a single layer in the Xi direction, there is only one coefficient a and therefore a single complex valued 
suplementary variable. The equations for the magnetic and electric fields are 

dtEi - (V A H)i - (jEi + 0^14^1 = 0, dtHi + (V A E)i - aHi + a^W2 = 0, 

dtE2 - (V A H)2 + o£:2 - 0, dtH2 + (V A E)2 + aH2 = 0, 

dtEs - (V A H)3 + 0^3 = 0, dtHs + (V A E)^ + aHs = 0, 

dtWi + oM^i - £^1 = 0, dtW2 + oM^2 - i?i = 0, 

In 2d this is identical to the layers in [27] and equivalent to those in [2]. 
The principal symbol and lower terms are 

Reversing the computation shows that (y, M^) S Ker i?(T,^) if and only ii V — MU, W = 1/tU, and 
Li{t, t^p^, • - - , T-^CT3 )U = 0. The characteristic polynomial is therefore the same as for Berenger's layer. 
Thus, by Theorem O 

deti?(.,0 =r^-E(7!^- 
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Theorem 3.6. If aj{xj) € L°°(R) and vanish for \xj\ < Lj then the system (|3.29p for V := MU and 

Wj := Vj /{t + aj{xj)) is strongly well posed in L^(]R'^) and perfectly matched in the sense that for sources 
supported in R := Ilj]{\xj\ < Lj} the function U computed from V, W agrees inM.t x R with the solution 
of Maxwell's equation with corresponding sources. 

Proof. The proof of Theorem 13.41 applies with only minor modifications. □ 

Remark 3.3. Note the ease with which strong well posedness is established and the lack of regularity 
required of the functions aj. 

4. Analysis of layers with only one absorption by Fourier-Laplace transform 

There are cases where the energy method presented above does not prove well posedness. This is the 
case for the Berenger algorithm when the ellipticity assumption is not satisfied and the absorptions are 
not regular. Notably for the Maxwell system and discontinuous absorptions. In this section we present a 
systematic analysis by Fourier-Laplace transformation of transmission problems with absorption in only 
one direction. 



4.1. Fourier analysis of piecewise constant coefficient transmission problems 

Return to the situation of (|1.3p with operators L and R on the left and right half spaces and transmission 
condition (|1.4p . Suppose that both L and R are weakly hyperbolic in the sense of Garding. An example is 
the classical method of Berenger with one absorption. Among other things we will prove that the method 
is well posed and perfect. Note the open problem at the end of the introduction emphasizing that we 
do not know if the classic algorithm with two discontinuous absorptions is well posed. In addition, we 
show, by a non trivial analytic continuation argument in ^4.1.5i that the perfection of Berenger's method 
can be verified using the modified plane wave solutions from his original paper. It is our hope that the 
analysis may help in the construction of new perfectly matched layers. 

4.1.1. Hersh's condition for transmission problems 

This section takes up the analysis of mixed problems following Hersh in |llj . In the present context we 
treat transmission problems which are essentially equivalent. The analysis of Hersh supposed the interface 
is noncharacteristic which is never the case for Maxwell's equations. We address the changes that are 
needed to treat problems with characteristic interfaces. 

First analyse the solution of the constant coefficient pure initial value problem £U = F on R^+'^ by 
Laplace transform in time and Fourier transform in x' ~ {x2, ■ ■ ■ ,Xd)- The transform 

U{T,xi,y) := j e~^* {27T)-'^/^e-'''='-''^ U{t,x') dtdx' 

decays as |xi| — )■ oo and satisfies 

C{T,d/dxi,ir])U = F in R. 

When Ai is invcrtible this is a standard ordinary differential equation in xi. When Ai is singular, the 
analysis requires care. The homogeneous equation C{T,d/dxi,iri) U = has purely exponential solutions 
gpxi corresponding to the roots p of the equation 

det£(T, p, irf) = 0. (4.1) 

Hyperbolicity of C guarantees that for Rer > tq and 77 G W^^^ this equation has no purely imaginary 
roots. 

The number of boundary conditions at xi = for the boundary value problem in the right half space 
is chosen equal to the number of roots with negative real part (see also Remark |4.1|) . That integer must 
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be independent of t, 77. Since roots cannot cross the imaginary axis, the only way the integer can change 
is if roots escape to infinity. That can happen when the coefficient of the highest power of p vanishes. 
The next hypothesis rules that out. 

Definition 4.1. A hyperbolic operator C{dt,dx) is nondegenerate with respect to xi when there 
is a Ti > so that the degree in p of the polynomial det C{t, p,iri) is independent of (r, 77) for Rer > 

Ti , 7] e K'*-i. 

Example 4.1. 

1. In the noncharacteristic case, det Ai 7^ 0, the condition is satisfied and the degree with respect to p is 
equal to N. 

2. For Maxwell's equations written in the real 6x6 form, the degree with respect to p is equal to 4. If 
written in the complex form p.2p . the degree is 2. 

3. The formula for the characteristic polynomial in Theorem 12.21 shows that if L is nondegenerate then 
so is the Berenger doubled operator L with one absorption cti in xi > 0. The degree in p is the same for 
L and L. 

4. If £ = Ci+B is nondegenerate with respect to xi then so is the operator P := £i{d)+a^^ B = a~^C{ad) 
for any a > 0. If the degree for £ is constant in Rer > ti, then the degree for P is constant for 
Rcr > a^^Ti. 

For the lemmas to follow it is useful to transform so that Ai has block form. 
Lemma 4.1. If C in (|2.ip is nondegenerate with respect to xi then for Rer > ti , £ M'^^^, 
(i) the degree in p of the polynomial dot £(r, p, irj) is equal to rank Ai , 

(ii) the number of roots p with positive real part is equal to the number of negative eigenvalues of Ai- 
Proof. 

Since C is nondegenerate, it suffices to study the case 7; = 0. 
(7) Choose invertible K so that 

K^^ Ai K = ( ^ ) , A an invertible square matrix of size rank^i . 







K ^C(t,p,0)K = (^^^"^^ + matrix independent of T, p . 

V T I / 



Then 



It follows that the degree in p is no larger than rank Ai . 

The coefficient of ^''^-ik.Ai ^^g-j- jg ^ polynomial in r of degree < N — rank Ai- For large t 

the coefficient is equal to 

( dot A) (r^-^^""^ ) + lower order in r . 
Thus the degree in p is rank Ai for such r proving the result. 



dct/:(r,p,0)=dct( [^^y-^ ^r]+B 



(ii) 

Hot r.(T n a) = Hot ( r ^ 

_ r/y 

For fixed r, sufficiently large, p 7^ 0, and p/r is a root of the polynomial p{x, l/r) of degree rankyl: 

p{x e) = ^ + + "^11 "^12 

£i?21 I + £^22 

This polynomial has exactly rank^ roots. By Rouche's theorem. 

Pi 1 

T Xj 
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where the Xj are the rank A eigenvalues of A repeated according to their algebraic multiphcity. Since the 
eigenvalues of A and the nonzero eigenvalues of Ai are the same, this completes the proof. □ 

Remark 4.1. For the transformed one dimensional hyperbolic operator C{dt,di,iri) the number of in- 
coming characteristics at the boundary xi = in the right half space is equal to the number of strictly 
positive eigenvalues of Ai . The second part of the Lemma shows that this is equal to the number of roots 
with negative real part. The two natural ways to compute the number of necessary boundary conditions 
yield the same answer. 

The next lemma shows that for nondegenerate operators, the characteristic case can be transformed 
to a standard ordinary differential equation. 

Lemma 4.2. Suppose that A, M G Hom (C^) and the equation det(^p + M) = has degree in p equal 
to rank A and no purely imaginary roots. Then, 

(i) The matrix M is invertible and all solutions of the homogeneous equation 

A^ + MU ^ (4.2) 

ax I 

take values in the space G := (Range A) satisfying dim G = rank A. 

(ii) There is a M g HomG so that a function U satisfies (j4.2|) if and only if U is G valued and satisfies 

dU ~ , , 

— + MU = 0. (4.3 
dxi 

(Hi) The vector space U of solutions of (|4.2p is a linear subspace o/C°°(K) with dimension equal to 
rank A. The Cauchy problem with data in G is well posed. 



Proof. 

{i) Since p = is not a root, M is invertible. The equation U = M^^AdU/dxi shows that continuously 
differentiable solutions U takes values in G. More generally, if J7 is a distribution solution and i/' G C'o°(I^) 
takes values in the annihilator, G^ of G, then 

(C/,i/') = {M-^ AdU/dxi, ^jj) = {dU/dxi, {M-'^ A)*^j) 

But G = range Af-M so G-^ = kcr(M-M)*. Therefore {M'^ A)*-4' = so (t/,V') = which is the 
desired conclusion. 



(a) Multiplying the equation by an invertible P and making the change of variable U = KV transforms 
the equation to the equivalent form 

PAK^ + PMKV = 0. 
dxi 



Choose invertible P, K so that PAK has block form 

PAK = 



0/ ' 



where / is the rank A x rank A identity matrix. With V = {Vi, V2) one has the block forms 



One has 

H2I H22 



dct{Ap + M) = detp-i detf^^ + ^"M detA^i. 
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The first part of tlie preceding lemma implies that the determinant on the left is a polynomical of degree 
rank A in p. It follows that H22 is invertible. 

The solutions V satisfy H21V1 + H22V2 = so take values in V :— {V2 — — iJ22^iJ2iVi}. The function 
is a solution if and only if it takes values in V and 

dVi _i 

— — + RVi =0, R :~ Hii — Hi2Hy2 H21 ■ 

ax I 

If : V V is the map, 

(Vi,V2) ^ (rVi, -HZ^H2iRVi), 



then y is a solution if and only if it is V valued and satisfies dV/dxi = NV . Writing V = K and 
M = -KN implies {ii). 

{iii) Follows from (ii). □ 

Lemma 4.3. If C is hyperbolic and nondegenerate with respect to xi, then its principal part Ci{dt,dx) 
is also nondegenerate with respect to Xi. The degree in p of dot Ci{t, p,ir]) is constant for Rcr > and 

Proof. With notation from the preceding proof, 

A.„*.,.,,^.-'((t-:)-g;!tt!S!';t!))-'- 

The proof of the last lemma showed that for Rcr > ti and rj e M''^^, H22{T,ri) is invertible. 

The computation in Lemma 14.11 shows that for t] — and R 9 t 00 the coefficient of p'^^'^^ -^1 
has modulus > cr^"''''"'' -^^ with c > 0. This implies that 77 = is noncharactcristic for H22- Therefore 
H22{dt,dx') is hyperbolic. 

Replacing C by its principal part Ci has the effect of replacing each operator Hij (d) by its principal 
part. This yields identity (|4.4p with C and the Hij replaced by their principal parts. 

Since the principal part of a hyperbolic operator is hyperbolic, it follows that {H22)i{dt,dx') is a 
homogeneous hyperbolic operator. Therefore {H22)i{T,irj) is invertible for rj E M''^^ and Rcr ^ 0. Thus, 
the coefficient of p'^'^^ -^i in dct £i(t, p, irf) is nonzero for 7] e W^-^ and Re t 7^ 0. □ 



Lemma 4.4. Suppose that the ordinary differential eguation (|4.2p satisfies the hypotheses of Lemma \Jf..S\ 
Denote by the linear space of solutions which tend exponentially to zero as Xi — > ±00 and by 
their traces at xi =0. Then 

(i) E^OkciA = {0}, 

(ii) d\m AE"^ = dim 

(iii) The map U i— > C/(0) is an isomorphism from E^ to E^ , 

(iv) AE+ n AE- = {0}, 

(v) AE+ ® AE^ = Range A. 

Proof. 

(i) The absence of purely imaginary roots shows that every solution is uniquely the sum of two solutions. 
One grows exponentially at +00 and decays exponentially at xi — —00. The second grows at —00 and 
decays at +00. In particular the only bounded solution is the zero solution. 

If e+ € E^ n ker A, denote by U{xi) the solution with this Cauchy data. The function that is equal 
to U on xi > and equal to in xi < is a distribution solution of ()4.2|) on all of M since A[J7]a;j=o = 0. 
This solution is bounded hence identically equal to zero. Therefore e+ = 0. The case for E" D ker A is 
analogous. 
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{ii) Follows from (i). 



(iii) It is surjective by definition. If it were not injective for E'^ , there would be a nontrivial solution 
U{x) exponentially decaying as xi — >■ +oo with U{0) ~ violating (i). 

{iv) The set AE~^ consists of the values AU-\.{Q) with U+ satisfying (|4.2|) and exponentially decreasing in 
xi > 0. If the intersection were nontrivial there would be a solutions U- decaying as xi — — oo so that 
AU^{Q) = AU^{0). The function V equal to U+ in xi > and U- in xi < is then a distribution solution 
for all xi exponentially decaying in both directions. Hyperbolicity implies that V ~ contradicting the 
nontriviality. 

{v) Using (m) and (iv), one sees that the direct sum on the left is a subspace of Range A of full dimensioD 

The next lemma is needed in in §4.1.2. 

Lemma 4.5. Assume that the hypotheses and notations of Lemma ^'^s force. Then for K G 
there is an F ^ C'o°(] ^ oo, 0[) so that the unique solution of 

A^ + MU ^ F, lim \\U{xi)\\ = 0, (4.5) 

satifies U{0) = K. 

Proof. Consider first the case of invertible A. A change of dependent variable yields the block form for 
the new variable still denoted U 

iki + (^0^ ^ = ^ = ('^i' ^2)' F = (^1' ^2) , specM± c {±Rez > 0}. 

Then E+ = {U2 = 0} so K = (i^i, 0). Choose F = (i^i, 0). Then [/(O) = K if and only if, 

Ki = / e^'+'Fiis) ds. 

J —OO 

This is achieved with, 

F,{s) = x{s)e-^'+'K,, ;^gCo°°(]-oo,0[), J x{s) ds ^ 1 . 

When A is not invertible change variable as in Lemma 14.21 to find the block form 

(lO\dU (HnHu\ ^ 
\0 Oj dxi ^ \H2i H22) 

with invertible H22. 

Part (i) of (|4.2|) implies that the map 3 G — {Gi,G2) 1^ Gi is an isomorphism. Write G 3 K = 
{Ki,K2). Choose F = {Fi, 0). Then choose a G valued solution U defined by 

^ + HnUi = Fi, U2 = -H^,^H2iUi. 
dxi 

One has U{0) = if and only if Ui{Q) ~ Ki. The construction in the invertible case finishes the proof. 

□ 



Suppose that 

L = dt + Aidi + --- and i? = + ^i^i + • • • 

are nondegenerate with respect to xi. For Rer > tq and 77 e M''"^, define E^ir^r/) to be the set of 
solutions of 

L{T,d/dxi,ii])V = 
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tending to zero as xi — > ±oo. Denote by E^{t,7]) C the Hnear space of traces at xi = of solutions 
in£;±(r,ry). 

Similarly with a possibly larger value still called tq, there are _E^(t, ry) C C^^ so that the solutions 
of R{T,d/dxi,ir])Z = taking values in E^{T,r]) are exactly those tending to zero exponentially as 
xi — ±oo. The subspaces Ef{T, rj) and E^{t, r/) depend smoothly on t, rj for Rer > tq and r/ G K''^^. 

Lemma |4. II implies that 

dim Ej^{t,i]) = # positive eigenvalues of Ai, dim E^{T,rj) = # negative eigenvalues of . (4.6) 
Consider the inhomogeneous transmission problem, 

LV ^0 when xi < 0, RW = when xi > , (4.7) 



{V,W)~g e Af when xi = . (4.8) 

The problem with inhomogeneous term F can be reduced to this form by subtracting on the left 
a solution of the hyperbolic Cauehy problem LU = F on R^+'' with C/|t<o = 0. Denote by 
V{t, xi,ri),W{T, xi,r])^ g{T,r]) the Fourier-Laplace transforms. The transform U is defined for xi e M, 
while V (resp. W) is defined for xi < (resp. xi > 0). The transforms V,W decay as |a;i| — > oo. V,W 
satisfy the ordinary differential transmission problem 

L{T,d/dxi,u])V = in xi < 0, R{T,d/dxi,iri)W = in xi > 0, (4.9) 

{ViT,0,r,),WiT,0,r,)) - 5(5,7?) e Af . (4.10) 

Hersh's necessary and sufficient condition for well posedness of the transmission problem is derived 
as follows. Uniqueness of solutions of (|4.9p . (|4.10p for Rer > tq, i] G R*^^^ is equivalent to the fact that 
there are no exponentially decaying solutions of the homogeneous transmission problem. That is, 

M n {EZ{T,v)^E+{r,rj)) = {0}. (4.11) 

In order to guarantee existence, one imposes the maximality condition, 

®{El{T,r])xE+{T,-n)) ==C^xC^^. (4.12) 

Using (|4.6p , this determines the dimension of Af from the coefficients Ai and Ai of L and R respectively. 

Definition 4.2. If the transmission problem ()4.7|) . (|4.8|) satisfies (|4.12|) for aU Rer > tq and r; e K''"^ 
it is said to satisfy Hersh's condition. 

Theorem 4.1. Hersh's condition is satisfied if and only if there is an r and a Aq so that for all A > Aq 
and g supported in t > with e^^*^g E H''^^{Rf^,) with values in x C^^ there is a unique V,W 
supported in t > with 

e-^V e i?"(] -cx3,cx3[x{xi < 0}) and e"^*!^ £ iJ^(] - 00, c»[x{xi > 0}) 

satisfying the transmission problem (|4.7|) . (|4.8p . 

Sketch of Proof. We have shown that the Hersh condition permits one to compute a candidate Fourier- 
Laplace transform. We outline how the condition implies the desired estimate. The method is to use the 
Seidenberg-Tarski Theorem 12.31 to derive a lower bound on the real parts of the roots v together with 
a contour integral representation. The same elements form the heart of [18], and §12.9 of [13]. In the 
present context we treat a transmission problem rather than a boundary value problem. In addition, one 
needs to use the earlier lemmas to treat the case when xi = is characteristic. 
Choose A > max{To(i), To(i?)}. The equations 



det L{T,i>,iri) = 0, dcXRir^v^ir]) = 
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with RcT > A, rj E M.'^"^ have no purely niiagmary roots. Define 

C(i?) := min||Rez/| : j^eW^'^ Rer>A, |rp + jryp < i?^ 

{ dot L(r, I/, = or det i?(r, i^,ir/) = 0}|. 
The Seidenberg-Tarski Theorem 12.31 imphes that there is a p € Q and 6 7^ so that 

C{R) = Wib + oiR)), as R^oo. 
Thus, there are C,N so that, for any r, 77 with Rer > A, 

|Rei.| > f (4.13) 
1 + |(r, ?7)P^ 

The solutions in (r, 77) are written using a contour integral representation of W in the block form 
of Lemma 14.21 Here the matrix Hij depend of (t, r^). Denote by D = D{T,ri) the finite union of squares 
with centers at the roots with Jicv < 0. The side of each square is the smaller of 1 and half the distance 
of the root to the imaginary axis. Then 

Wi = ^ i e^^i (r + {Hii + Hi2H22^H2i)) ' dr Wi, Wu = H^2^H2i Wi . (4.14) 



27ri JQD 

The Seidenberg-Tarski Theorem 12.31 applied to 

max||i(;|2 : \zf = l, H22W ^ z, Rcr>A, \Tf + \if<R'^^ 

proves that 

\\H22{T,r])-'\\ = RP{a + o{l)), a ^ 0, /3 G Q . 
This estimate together with (j4.13p yields with new C,N, 

\WiT,x,,r,)\^ dx, < C(l + |(r,r;)|2^)|W?,(0)|2. 

With the analogous expression for V the solution of (|4.7p satisfies 

V{T,xuv)f dx^ < C{l + \{T,f^)D\Vim\ 

The Hersh condition asserts that for each (r, 77), W7(0) and Vi{0) are uniquely determined by g{T, rj). 
Seidenberg-Tarski Theorem 12.31 yields an estimate 

\\Wi{0) , Vi{0)\\ < Cil + \{T,rj)\r\\g{T,rm'- 

The last three estimates together with Parseval's identity proves the desired estimate, 

3C,N, V.g, VA>A, \\e~^'U\\l < C 



L2(Ki + (d-i)) ■ 
\a\<N 



This estimate proves the existence part of the Theorem. □ 



4.1.2. Necessary and sufficient condition for perfection 

The Fourier-Laplace method is used to derive a necessary and sufficient condition for perfection of an 
absorbing layer. Begin with a closer analysis of the transform, J7(t, xi, 77), of the solution of the basic 
equation (|1.1[) . 

When Ai is invertible, U is analysed as follows. Denote by n±(T, 77) the projectors associated with 
the direct sum decomposition E^{T,ri) E£{T,ri) = C^. Define S±{t, xi,ri) as the Hom(C^) valued 
solutions of 



L(r, d/dxi, Z77)5± = 0, S±\^^^^ = A^'n± 
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Then S± decays exponentially as xi — > ±oo and 

X]-oofilS- + X[0,oo['5'+ 

is the unique tempered fundamental solution of L{t, d/dxi, irj). Decompose F = F- + _F-f, U = U+ + U- 
according to E^{t,vi) ® E^ir^ff) = . Then [/_ is the convolution of F_ with X]-oo.o[S- and [/+ is 
the convolution of F+ with X[o,oo['S'+- In particular t/_(r, 0,7?) vanishes on a neighborhood of [0,oo[ so 
C/(t, 0, 77) = U+{t, 0, rj) € E^{t, r]). The value of ?7 in xi > satisfy the homogeneous ordinary differential 
equation L{T,d/dxi,i7])U = 0. with initial value U{0) S E^{T,ri). 

To reach the same conclusion when Ai is singular, apply the lemmas of the preceding section to the 
equation L{T,d/ dxi,irf)Z = 0. Lemma [4.41 applied to A = Ai and M = r/ + i^j^2^jVj shows that 
both E^{t, 77) are subspaces of G and that the space of solutions is a direct sum E^{t, tj) © E~[{t, if). It 
follows that 

El{T,r])®E+{T,ri) = G(r, 7y) . 
Repeating the analysis in the nonsingular case applied to ()4.3|) shows that U{t, 0,7]) e E^{T,ri). 

Definition 4.3. For a transmission problem {L,R,J\f) satisfying Hersh's condition (Definition I4.2[) . 
Rer > To and 7/ £ R'^"-'^, define the reflection operator, H(T,rj) : E^{T,r]) —J- E'[{T,r]) as follows. 
Hersh's condition implies that for each K E E^{T,ri) there is a unique {V,W) S E2{T,r]) x E'j^{T,r]) so 
that {K,0) = {V,W)TtioAN. Define, H{T,r])K := V . 

Theorem 4.2. Suppose that the transmission problem (L,R,Af) satisfies the Hersh condition. The fol- 
lowing are equivalent. 

(i) The transmission problem is perfectly matched in the sense of Definition \l.l\ 
(a) There is a tq G R so that for all Rer > tq and r] € R''^^, F[{t, rj) = 0. 
(Hi) There is a tq E R so that for all Rer > tq and r] E R''^^, 

e E+{t, 7]), 31 Kr e E+{t, 7]), such that {Kl ,Kr)eN. (4.15) 

Proof. Conditions iii) and {Hi) are clearly equivalent. 

For the equivalence with [i), compare the values of U and V in {xi < 0}. Since both satisfy LZ ~ F 
and decay as xi — > — oo it follows that L{V — U) ~ 0, so, V — U := F is an Ej^ valued solution of LT ~ 0. 
Since F = in xi > 0, W E E^. The transmission condition requires that 

Af 3 {V{0),W{0)) = (&(0) + F(0),t?(0)) = (&(0),0) + (F(0) , t?(0)) . (4.16) 

Since (F(0), 11^(0)) E E~[{t, 77) x E^{t, 77), (|4.16p expresses {U{0), 0) as a sum of an element in M and an 
element of EJ^ (''': ^7) ^ (''': ^7) ■ The Hersh condition (|4.12p asserts that such a decomposition is unique. 
Therefore {V{0),W{0)) is uniquely determined from U{0). 

The method is perfectly matched if and only if for all F supported in xi < 0, i > 

V ^ U\ 

This occurs if and only F vanishes for xi < which holds if and only if F(0) — 0. 

If the method is perfectly matched, then in the decomposition (|4.16|) one has F(0) = 0. Then 
{U{0),W{0)) E M. Lemma 1431 asserts that for any K E E^ tl^cre is an F so that U{0) = K. This 
proves that (|4.15|) holds. 

Conversely if (|4.15p holds, then in the decomposition (|4.16|) . F(0) = so F = 0. It follows that 

Remark 4.2. 1. When (|415)) holds, the decomposition of {K, 0) £ x C^^ in the direct sum (|4T2)) is, 

{K,0) = {K,W(K)) - {0,W{K)) E Af®{E^xE+). 
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2. With K = U{0) as above the solution (V, W) of the ordinary differential equation transmission problem 
is given hy V = U\^<_o and W is the solution oi RZ = with Z{0) = -W{K). 

3. In the important case where TV = M, invertible Ai and Ai and transmission condition Af — {V = W}, 
the perfection criterion (Hi) asserts that E^{t,7]) = E^[t,vi). 

We present a typical example showing that the natural absorbing layers are virtually never perfectly 
matched in dimension d > 2. 

Proposition 4.1. Consider the dissipative symmetric hyperbolic example with d = N = M = 2, 
^i=Q_°^), A2=(^JJ), R = L + P, P = P*>0, JV^{iV,W):V = W}. 

(i) The transmission problem is perfectly matched if and only if P ^ 0. 

(ii) The corresponding problem with d — \ is perfectly matched if and only if P is diagonal. 



Proof. Define 



Ml := A-^^[t + zr/Az] , Mr A^^ir + ^77^2] + A^^P , 



so that A^^ L(t, di,iri) — di + ML{T,ri) and similarly for A^^R{T,di,iri). For Rer > and 77 £ K, the 
matrices Ml and Mr have one eigenvalue with positive real part and one with negative real part. The 
eigenspace corresponding to positive (rcsp. negative) real part eigenvectors is equal to ij^(r, 77) (resp. 
E~l^{T,ri)). Therefore the necessary and sufficient condition for perfection is that for Rer > tq and any rj, 
E+iT,zr,) = E+{T,ir^). 
Since 

, . I T + p ill \ 1 T I \ ' 2 2 

L(T,p,iri)~\ . , and, det L(t, p,ir]) ~ t" — p +77, 



the eigenvalue of AlL{T,ri) with positive real part is p = y'r^ + rp . The eigenspace is the kernel of 
L{t, p,iri). Therefore 

Et{r,v) - C{-tV,r + p). (4.17) 

Since Mr = Ml + A]\^P, a necessary condition is that the family of vectors v{ri,T) {—iij,T + p) 
be eigenvectors of the constant matrix A'^^P, which is possible only if A^^P is a constant multiple of 
the identity. Therefore P = cA\. Since P > and A\ has eigenvalues of both signs, it follows that c = 
proving (i). 

In the one dimensional case there is just one eigenvector (0, 1) which must be an eigenvector of A'^^P . 
Since (0, 1) is also an eigenvector of A^ it follows that (0, 1) must be an eigenvector of P . Since P = P*, 
the orthogonal vector (1,0) is also an eigenvector and P is diagonal. Conversely, if P is diagonal the 
condition is satisfied. □ 



Remark 4.3. 

1. Examples verifying perfection for a family of absorbing layers related to but not including those of 
Berengcr arc presented in [3]. To our knowlege Hersh's criterion for Bcrcngcr's layers has not been verified 
before. 

2. The perfection criterion is related to the plane wave criterion of Berengcr. We examine the relation in 
§4.1.6. 



4.1.3. Hersh's condition for Berenger's PML with piecewise constant ai 

Of our earlier results, only those of Section [3^ apply to discontinuous absorptions. So, if the generator is 
not elliptic, (for example, the PML Maxwell system of Berengcr), the preceding results do not prove that 
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the initial value problem is well posed. In this section we prove that the doubled operators of Berenger 
define a (weakly) well posed initial value problem provided that 

crj=0 for j>2 and cti (xi )= cr^ in , (4.18) 

and, the constant coefficient operators L on Rj. are both (weakly) hyperbolic. 

The unknown U satisfies (jl.Sp . Denote by = {U^, . . . , U^} the restriction of the unknown U to 
R^. They satisfy differential equations in the half spaces R^. 

Lemma 4.6. For U locally square integrahle on a neighborhood of {t,x) £ {xi — 0} the following are 
equivalent. 

(i) LU e on a neighborhood of{t,x) in Ri+'^ in the sense of distributions. 

(ii) There is a neighborhood O of {t,x) so that LU^ is square integrable on O C] R^ and [AiU\ = 0. 
Remark 4.4. 

1. The first hypothesis is often verified by combining LU + B{x)U G Lf^^, U G Lf^^ and B £ Lf^^. 

2. [AiU] makes sense since the differential equation implies 

9i(IiC7+) e iL(]0,e[; Hr^liRt,,)). 

With U G L2(]0,e[; i?,;;^(R'^) this implies that AiU+ G C([0,e[; H'^^iR'^ j). An analogous resuh holds 

— — —1/2 

for AiU~ . Therefore the traces from both sides and the jump are well defined elements of Hi^J . 

3. is clear on a formal level since if AiU were discontinuous there would be a 5(xi) term from the 
differential operator L applied to U . 

4. The standard proof based on these remarks is omitted. 

We have supposed that the nonzero data arc initial values f^{x). By the usual subtraction one can 
convert the problem to one with homogeneous initial values and right hand side and inhomogcneous 
transmission condition. In this way, the determination of is reduced to finding satisfying the 
inhomogcneous transmission problem 

Li{dud.,)W^ + B^W^ = 0, M\W] = g, (4.19) 

where 



/a^lN . . . 0\ 
... 



B^ 



(4.20) 



V ... 0/ 

and g{t, x) take values in Range Ai. The unknowns W and source g are vectors of length dN. 

Theorem 4.3. Suppose that L{d) is a hyperbolic operator nondegenerate with respect to xi and that the 
Berenger's doubled operator L is weakly hyperbolic for a ~ . Then, the transmission problem (|4.19p 
with absorption (|4.18p satisfies Hersh's condition. 

Proof. Drop the tildes on the Fourier-Laplace transforms of VK, g for case of reading. The transformed 
problem is 

(Ii^ + Z(t,0,z7?))w^ =0, M[W] = g. (4.21) 

The condition of Hcrsh is that for an arbitrary right hand side g in range Ai this transmission problem 
has one and only one solution. 



The uniqueness of solutions of (|4.2ip is equivalent to 



Denote by E~{T,r],a) the spaces associated to the Berenger operator operator L with absorption a. 



AiE^{T,ij,<j-) n AiE±{T,Tj,a+) = {0}. (4.22) 
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Existence is equivalent to 

AiE~{T,r),a-) + AiE±{T,i],a+) = rangeli. (4.23) 

Part (m) of Lemma 14. 41 implies that 

dim (r, r/,cr^)) + dim (yliij^(T, ry, a^)) = dim (range ^li) , 

so (14221) implies (|423l) . It remains to prove (|422|) . 

For (T > 0, the split Berenger operator L is hyperbolic so for Rer > tq{(7) the solutions of 

L{T,d/dxi,i7])U = (4.24) 

are generated by exponentially growing and exponentially decaying solutions. The next lemma identifies 
these solutions in terms of the corresponding solutions of 

L{T,d/dxi,iT])V = 0. (4.25) 

The result shows that the traces at xi = 0, E^, are independent of a. 

Lemma 4.7. For cr > 0, Rer > To{a), rj £ R''^^, 
(i) The map 

V{xi) ^ (^^AiT>((T + fT)xi/r) , ^A2i/((T + a)2;i/T) , ... , ^Adi/((T + a).Ti/r)^ 

is an isomorphism from solutions of (j4.25p onto the solutions of (|4.24[) . 
(a) fi is a root of det L{t, -Tirj) ^ if and only if v = {t + (T)/i/T is a root of det L{t, -^rf) = 0. 
(Hi) For the roots in {ii), the real parts of fi and v have the same sign. In particular, the map in (i) is 
an isomorphism Ef^ (r, rf) i— > E-^ (r, cr) . 

(iv) The map W = (Wi, . . . , Wd) ^ '^^ isomorphism E^(t, rj, a) — > E^{t, r/). 

Remark 4.5. In (i) it important to know that the solutions V{xi) are entire analytic functions of Xi 
so it makes sense to evaluate V at points off the xi-axis. In the literature this is sometimes called a 
complex change of variables. It is only reasonable for analytic solutions. A related idea is used in the 
Fourier-Laplace analysis for general cri(a;i) presented in §4.2. 

Proof of Lemma 14.71 

(i) U U = {Ui, . . . , Ud) satisfies (|4:24| then with W := Y.j Uj, 

+ {r + a)Ui ^ 0, tUj + irjjA^W = 0, j = 2,...,d. (4.26) 
Multiply the first by r and the last d — 1 by (t + cr). Sum and then divide by r to find, 

Ai^ + ^±^L{T,0,tTj)W^ = 0. (4.27) 
dxi T 

Conversely if W satisfies (|4.27p and Uj for j > 2 is defined from W using the last equations in (|4.26p and 
Ui:=W -J2j>2 Uj then U satisfies g^J). 

The solutions W to (|4.27p are exactly the V{{t + a)xi/T) with V satisfying (|4.25p . This proves that 
the mapping in (i) is surjective. 

The set of solutions V of (|4.25p has dimension rank Ai . The set of solutions of (j4.24p has dimension 
rank Ai = rank Ai (see ()2.3p ). so surjectivity implies injectivity. 



(ii) and (iv) follow from (i). 
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{Hi) Denote by K the mapping from (i). When t > 0, (t + <7)/t is also positive and real Therefore K 
maps decaying (resp. increasing) solutions to decaying (resp. increasing) solutions. Thus for t > tq and 
real, 

K{E+{T,^))=E±{T,r,,a). (4.28) 

For all Re r > tq, K{E^{t^ rf)) is a subspace of solutions of ()4.24|) with dimension equal to dim E^{t^ if). 
If (|4.28p were violated, K{E~^(T,r])) would contain exponentially growing solutions. If this happened at 
T, 77 with Re T > To consider the values r(r) = Re t + rIm t for < r < 1. For r = (|4.28p is satisfied 
while for r = 1 it is violated. Let 

Re — — fi : dot L{T{r) , ^ , i-q) = , Re /i < 0| . 

Then /(O) < 0, /(I) > and / is continuous, so there is a < r < 1 so that /(r) = 0. Then for 
T = Rer + r Imr there is a purely imaginary root. This violates the hyperbolicity of L establishing 
(|4.28|) This proves [Hi) completing the proof of the Lemma. □ 

We now finish the proof of Theorem 14.31 bv proving ()4.22|) . Lemma [4.71 implies that the spaces of 
Cauchy data E^ are independent of cr. Therefore if (|4.22p is violated then also 

Iii;~(T,r/,a+) n Ii^±(r,77,a+) ^ {0}. 

This contradicts part [iv) of Lemma 14.41 for the operator L with absorption cr+. The proof of Hersh's 
condition is complete. □ 



In these problems with only one nonzero absorption coefficient cri and cri = when xi < one can 
consider a transmission problem which is only split in xi > 0. The next result shows that this partially 
split problem satisfies Hersh's condition if and only if the fully split problem does. 

Introduce the partially split problem {L,R,M) where 

L^Li{d), R^Li + B+, with a+ > 0, JV := {(y,W) : V -^Wj € kerAi}, (4.29) 

j 

with 5+ given by (|4.20[) and the split variable on the right is W = {Wi, . . . , Wd). 

Corollary 4.1. Suppose that Uj — for j > 2, and > 0. Then the partially split Berenger transmission 
problem (Li, Li + ^N) defined by (|4.29p satisfies Hersh's condition if and only if the fully split problem 
does. 



Proof. Denote by {V, W) = {V, Wi, . . . , Wd) the variables for the partially split problem and [U, W) = 
(^{Ui, . . . ,Ud), {Wi, . . . ,Wd)) the split variables. If C/(t, xi, 77), PF(t, xi, 77) is an exponentially decaying 
solution of the split Laplace-Fourier transformed homogeneous transmission problem, then (1/, W) — 
Uj , W) is an exponentially decaying solution of the partially split homogeneous transmission problem. 
Conversely, if V(T,xi,ri),W{T,xi,r]) is a solution of the homogeneous partially split problem, the 
computation leading to (|4.26|) shows that 

Ui := -T-^AidiV, Uj := -r^^ irjjAj V, j > 2, 



is an exponentially decaying solution of the fully split homogeneous transmission problem. 

Therefore, if either problem has decaying solutions for rj real and Rer arbitrarily large, then so does 
the other. □ 
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AAA. Perfection for Berenger's PML with piecewise constant cti 



Theorem 4.4. With the hypotheses of Theorem \4-3\ the Berenger transmission problem is perfectly 
matched. The Berenger transmission problem that is only split on the right is also perfectly matched. 

Proof. Verify condition (m) of Theorem 14.21 For K e E~{T,ri,a^) consider the unique decomposition 
guaranteed by the Hersh's condition, 

{K,0) = iW-,W+) + {F-,F+), (4.30) 

where, 

{W-,W+) = (^{Wr,...,W^),{Wt,---,W+)) gAA, {F~,F+) e E^{T,v,a-)x E+{T,7j,a+). 

Perfection is equivalent to F^ = 0. 

By inspection, one such decomposition (|4.30p is given by 

(i^,0) = iK,K) + (0,-/0 

where we use the fact from Lemma 14.71 that 

E±{T,Tj,a-) = E±{T,r^,a+). 

As this decomposition satisfies F^ = the proof of the first assertion is complete. 

For the partiaUy spht case, K e E^{t, 77) has a unique decomposition from the Hersh's condition, 

{K,Q) = {W-,W+) + {F-,F+), 

with 

[W- ,W+)^{W~ ,{W+,..., W+)) e AA, {F- ,F+)eE^ (r, V) x E± (r, r?, a+) . 

Define 

Part {i) of Lemma Wl\ implies that e ''l^'^^)- addition, J2j = so {K, W^+) G TV. 

By inspection 

{K,0) = (A' , W+) + (0 , ~W+) 
is the unique Hersh decomposition. Since F^ vanishes for this one the proof is complete. □ 

4.1.5. Analytic continuation for Maxwell like systems and Berenger's plane waves 

In this section we investigate Berenger's method for operators, including the Maxwell system, whose 
characteristic polynomial is t'p{t'^ — I^P)''- For ease of exposition we treat the case d = 2 and the explicit 
operator. 

Analogous results are valid for the Maxwell system with only slightly more complicated formulas. 
For Re T > and r/ e M there is exactly one root of det £1 (t, p, irj) ~ with Re p > given by 

p = \/t^ +if, Rcp>0. (4.32) 

The corresponding eigenspace E'1{t,vi) from (|4.17p is spanned by <J>(t, 77) = {—iri,T + p). If Li is 
the Berenger operator doubled in the xi direction one has the same roots and E^ is spanned by 

(pAi$,iryA2$). 



Proposition 4.2. 
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(i) For each rj, pir^rj), EJ^[T^rj), and E'^{T^rj) are holomorphic in Rer > with continuous extension 
toRcT>0. 

(ii) If a I > then for Rer > and 7^ G R the equation det i(T, 177) = has exactly one root v with 

positive real part. It is given by v = (r + (Ti)p/t. 
(Hi) For ai > 0, the relation (j4.12p with Li on the left and L on the right is satisfied on {Re t > , p 7^ O} . 
(iv) The mapping H{T,ri) is for each rj holomorphic in Rer > with continuous extension to {Rer > 

0,P^0}. 

Proof, {i) For Rer > there are two roots ±p with p from (|4.32|) . One has strictly positive real part 
and the other strictly negative. Each is holomorphic in Rer > 0. Holomorphy for $(r, rj) follows from its 
expression in terms of p. As $ is a basis for E^{t,7]) holomorphy of the latter follows. 

So long as the eigenvalues ±p remain apart as Re r — >■ they and their cigenspaccs will be holomorphic. 
The delicate case is when r^ + 77^ — 0. The limiting points are (±177, r/). 

If 77 = 0, then p = T and the eigenspace is (0, 1). Both arc continuous up to the boundary. 

When rj ^ one has p — >■ so p is continuous up to the boundary. Then $ is continuous up to 
the boundary and nonvanishing from its expression in terms of p. Therefore E'^ is continuous up to the 
boundary. 

(a) It suffices to remark that this is an eigenvalue and then to show that the real part is positive. For 
the latter compute 



da 

For Re r > this has positive real part so the real part of the eigenvalue is increasing as a function of a 
so is positive for all cr > 0. 



{Hi) It suffices to show that (|4.1ip is vahd for such 5,77. Suppose {v,w) = (ui, W2, wi, W2) <S x E^. 
Must show that vi + V2 7^ wi + W2. Since {wi,W2) G E^ it follows that wi + W2 € E'^. Similarly 
V1+V2 G E^-^ ■ Thus it suffices to show that i?^^ and E'j^ are uniformly transverse as Re r — > 0. It suffices 
to show that (177, t + p) and (177, t ^ p) are uniformly independent. This follows from p 7^ 0. 

{iv) The holomorphy of H follows from (i). The continuous extension follows from (i) and {Hi) . □ 

Since the method is perfectly matched, H ^ for Re r > 0. By continuity the map vanishes for purely 
imaginary r 7^ 0. This shows that for {Rer > 0} \ 0, the function equal to 

e'^*+P(^''iy-'^+'1'-'^^ for Xi < 0, and e'^t+Pir,r,)x,+^nx, ^-apx^jr ^ fQj. 

satisfies the Berenger transmission problem. For Re r > these solutions decay (resp. grow) exponentially 
as Xi — > 00 (resp. x\ -> —00). Though such solutions serve to verify perfection they don't look very 
physical in isolation. 

On the other hand, when r is purely imaginary and not equal to zero, the solution is a bounded plane 
wave in a;i < and is a plane wave modulated by an exponentially decaying factor in x\ > 0. These are 
the solutions which Berenger constructed to show that the method was perfectly matched. 

In the language of the analytic objects constructed in the preceding lemma, Berenger's plane wave 
solutions show that H{is, 77) = when s is real valued with > t]^. For 77 fixed the function r 1— >■ H{t, rj) is 
holomorphic in the right half plane continuous up to the imaginary axis punctured at ±i|77|, and vanishes 
on the boundary interval r = is G iM with > 77^. By Schwarz reflection and analytic continuation this 
implies that H vanishes in the right half plane. 

In summary, the computation of Berenger is actually sufficient to prove perfection for Maxwell 's system 
given the structures provided in this paper. 

Remark 4.6. The perfection argument based on plane waves is not valid in full generality where the 
objects like E and H are analytic in Re r > ro with tq > 0. This is the case, for example, whenever the 
absorbing layer is amplifying. 
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4.2. Fourier- Laplace analysis with variable (Ti(xi) 

Consider the case of only one nonzero CTi (xi ) . If L is hyperbolic for one constant value 7^ the scahng 
{t, x) I— >■ {at, ax) shows that L is hyperbolic for ai/a. Therefore L is hyperbolic for all constant values cti. 

The results of ij4.1l will be extended to the case (Tj = for j > 2 and variable coefficient cri(xi). The 
Fourier-Laplace transform C/(r, xi,?;) of the Bcrengcr split operator satisfies 

L{t, d/dxi,r])U = F, — 00 < xi < 00 , 

with variable coefficient cri(a;i). 

The first line of the proof of Lemma 14771 vields ()4.26|) with a = ai{xi). As in the proof of that lemma 
one derives ()4.27|) now with a = ai{xi). The important observation is that the xi dependence of the 
coefficient appears only as a scalar prefactor in (|4.27|) . Such equations will be analysed in the same way 
as the equations in Lemma 14.21 



4.2.1. Well posedness by Fourier- Laplace with variable <7i{xi) 

Theorem 4.5. Suppose cjj = for j > 2 and cri(a;i) G is real valued. Suppose in addition 

that L is nondegenerate with restect to Xi, and for one value ^ 0, L is hyperbolic. Then there is a 
To > and m so that for all A > tq and F S e'^*L'^(M.t '■ ^'"(J^t'a;')) there is a unique solution solution 
U £ e'^* L'^ {M.'^'^^) to the Berenger split problem LU — F . In addition, there is a constant C independent 
of F, A .so that 

< C^|h-"i^lL.(,^.,.(«.^,))- (4.33) 

Remark 4.7. 

1. The condition U € e^*L^ implies that U tends to zero at t — > —00 as does F. 

2. If F is supported in t > to it follows from (|4.33|) on sending A — > 00 that U is supported in t > io- 

Proof. The values of the Fourier Laplace Transform oi W = J2 ^i'^ computed from the ordinary 
differential equation 

A,^ + ^±^X,(r,0,z,) W = F. (4.34) 
axi T 

As in Lemmas 14.21 and 14. 3[ transform to the equivalent form, 

Denote the decomposition as = (Wi, Wn) and similarly F. The invertibility of H22 from Lemma 
yields, 

Wn = H^2'{Fii -H21W1) . (4.35) 
It suffices to find Wj which is determined from, 

dW^ + + Mir, rj)Wi = G , A/(r, ry) := Hn - ^12^/22' ^21, G := Fj + H-^Fu . 

dxi T 

The hyperbolicity of L implies that M has no purely imaginary eigenvalues. Correspondingly there is 
the decomposition, into the spectral parts with positive and negative imaginary parts, 

Wi = W+ + Wj-, G = G+ + G-, M = M+ ® A/-. 

For a constant, part {Hi) of Lemma 14. 71 (using the hyperbolicity of L) implies that for Rer sufficiently 
large (depending on cr), one has the spectral decomposition, 

M{T,ri) = M(t,77)+ ® M{t,t])- 
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corresponding to spectra with positive and negative real parts. 

Lemma 4.8. If g{xi) satisfies dg{xi) / dxi = cri(a;i). Then 

\dxi /V / \dxi T J 

Proof of Lemma. Since {de3^^^''Ui)/dxi = e9'''l''{g'MUi/T + dUj/dxi) one has 

\dxi /V / \dx\ ^ T T ' J 

proving the desired identity. □ 

Therefore 

The unique fundamental solution of di + M is equal to, 

e""'^^^X[o,oc[(a;i) + e-^^^^' x]-oo.,o\{xi) . 

Therefore, 

e^^^/^WT = (e--*^" X]-oo,o](xi))*(e^^^/^G)-. 
The kernel of the integral operator mapping to Wj' is equal to, 

exp(^- (xi -yijj^ Af(r,?7)+JJ Xxi>j/i • (4.36) 

Lemma 4.9. 

3To=To(Ai), VReT>To, V77eR'^"\ Vae[-/^,/i], spec ^— t^A/+(r, r/) C {Rez>0}. 

T 

Proof of Lemma. Part (iii) of Lemma H771 allows one to choose t\ so that for a = [i one has the desired 
conclusion for Rer > t\. Then for A € spec Af(T, 77)+ one has 

ReA>0, Re(l + -)A = Re^-ti^ > 0. 

T T 

For < c < write cr = a + b^ with nonnegative a, b summing to 1. It follows that Re (1 + /i/T)A > 0. 
This proves that ri suffices to treat the nonnegative values < cr < ^J.■ 

Choosing T2 for a = — /z, that value suffices for — ^ < cr < 0. Set tq equal to the maximum of ti and 

T2. □ 



The Seidenbcrg-Tarski Theorem 12 .31 shows that the absolute values of the real parts of the eigenvalues 
of A/(t, r/) are bounded below by C{\t\ + for some A^. And also that the spectral decomposition 

V I— >■ {V'^,V~) and its inverse are both bounded polynomially in |t, ryj. More generally for T,ri,fi,a as 
above, 



T - 

spec 



±^M+{t,v) C {|Rez| >G(|r| + H)-^}. 

Taking /i := one finds that for all xi,yi, the matrix 

r+{9{x,)-giy,)/ix,-y,) ^^^^^ ^^^^^^ 

T 

has spectrum in 

{Rez>G(|r| + |?7|)-^}, C>0. 
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The elements of the matrix (|4.37|) are bounded above polynomially in \t,t]\. Therefore the kernel (|4.36|) 
is bounded above by 

|r, 77|^ cxp ( - c(xi - yi)/|T, 77!^) x.^>y^ , c > . (4.38) 

This is proved using Schur's Theorem to reduce to upper triangular form by orthogonal transforma- 
tions of the spectral subspaces. Then solve the differential equation X' + M^X = by back substitution 
to prove |lexp(pM+)|i < C|r, ?7|Pe~'=''/l^'''l" . 

The operator with kernel (|4.38p is convolution by an element of L^(R) whose norm grows poly- 
nomially in |t, Tyj. By Young's theorem one concludes that the operator with kernel (|4.36[) has norm in 
Hom(L^(]R)) which grows at most polynomially in |r, 77I. 

There is an entirely analogous estimate for the expression for the spectrum with negative real part. 

Therefore, 

||M?Kt,xi,77)|L.(r) < Ci(H-|r|-f |r;|)^||G(r,xi,,7)||^,(^^ < C2 (1 + Irj-f |r;|)^ ||F(t,xi,,7)|L.(r) ■ 

A similar estimate for Wu follows from (|4.35p . Estimates for Uj follow from the second equation in (|4.26p . 
Plancherel's Theorem then implies (|4.33p . proving the existence part of well posedness. 

Uniqueness is proved by a duality argument of Holmgren type using existence (backward in time) for 
the adjoint differential operator (details omitted). □ 



4.2.2. Perfection for Berenger's PML with variable coefficient ai{xi) 



Lemma 4.10. Suppose that A, M satisfy the hypotheses of Lemma \4.S\ with G and M G Hom G are from 
that Lemma. Suppose in addition that f e L^^(R; C) and g is the unique solution of 



0, 



fis) ds. 



Then for 7 G G the unique solution of the equivalent initial value problems for the G valued function U , 



+ f{xi)M [/ = 0, cquivalently, 4- f{xi) M U 

dxi dxi 



0, 



C/(0)=7, 



Uix,) 



Proof. Compute using the differential equation, 



d 
dxi 

The lemma follows 



g(xi)M jj 



\dxiJ 



M ■ 



dU 
dxi 



fM^fM 



□ 



The next result shows that when the Berengcr split problem with absorption a{xi) defines a stable 
time evolution, then the problem is perfectly matched. Either the split problem is ill posed, or it is well 
posed and perfect. 



Theorem 4.6. Suppose that cri^x) £ L° 



has support in [0, p] for some p > 0, that aj = fo 



j ^ 1, and that the operator L with these absorptions is nondegenerate with respect to xi and defines 
a weakly well posed time evolution. Then, the L evolution is perfectly matched in the sense that for 
F S CQ°({t > 0} n {xi < 0}) the solutions U and U' with and without absorptions respectively, 



LiU 



F. 



U\ 



0, 



and, 



satisfy 



U 



xi<0 



U'\ 



LU' 



xi<0 ■ 



F, 



U'\ 
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Proof. Denote by U and U' the Fourier-Laplace transforms. The functions are characterized by 

ii(r, d/dxi, ?y) [/ = F, and, L{T,d/dxi,ri)U' ~ 

both required to decay exponentiahy as — >■ oo. The strategy is to construct a solution of the problem 
defining U' from the solution U. 

The equations for W = Uj and W ~ U'j in xi> have the form 

, dW — - , dW' T + ai(xi) ^, 

Ai— + MW = 0, Ai- — + ^-^MW = 0. 

dxi axi T 

Lemma [4. 101 applies with f{x) := (r + a{xi))/T. 

Define g as in that lemma. Set V ~W va xi <Q. For 2^1 > define 

The resulting function satisfies the differential equation required of W. In addition since e"^*^^^-'*^ is 
independent of xi for xi > p, V decays as rapidly as W. Therefore V satisfies the conditions uniquely 
determining W' . Therefore V = W' , and W'\^,<o = W^Ui<o- Use (jTM]) to recover U , U' from W , W' 
shows that U'\xi<o = t/|2-i<o, proving perfection. □ 



Example 4.2. 

1. If Li{Q^dx) is elliptic then CoroUarv 13.21 shows that the evolution of L is strongly well posed. This 
includes the case of anisotropic wave equations for which the layer is amplifying showing that perfection 
is not at all inconsistent with amplification. 



2. For the Maxwell equations and cri(a;i) € well poscdness is proved in the remark following 

Theorem 13.31 and we deduce perfection. 



5. Plane waves, geometric optics, and amplifying layers 

This section includes a series of ideas all related to plane waves and short wavelength asymptotic solutions 
of WKB type. We first recall the derivation of such solutions from exact plane wave solutions by Fourier 
synthesis. Then we review the construction of short wavelength asymptotic expansions. These are then 
applied to examine the proposed absorption by the aj . In many common cases the supposedly absorbing 
layers lead to asymptotic solutions which grow in the layer. Related phenomena arc studied by Hu, and 
Becache, Fauqueux, Joly [14], [5]. For the Maxwell equations for which the PML were designed, the 
layers arc not amplifying. At the end of Section 15.31 situations where the amplification docs not occur 
are identified. 



5.1. Geometric optics by Fourier synthesis 

When the coefficient a vanishes identically, both L and L are homogeneous constant coefficient systems. 
When (r, ^) is a smooth point of the characteristic variety, denote by r = r(^) the smooth parameteriza- 
tion, and I{l{t,£) and I{j^{t,£) the associated spectral projections for ^ ~ ^, see p.8|) . The function r(^) 
is homogeneous of degree 1, while the projectors are homogeneous of degree 0. The next argument works 
equally well for L and L. 

For G(^) G C^{W^) construct exact solutions for < e << 1, 

U'{t,x) := j e^^--+<^^'^UL{T,i)G{i-i/e)d^. 

Make the change of variable 
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and extract the rapidly oscillating term e'^^i'^^-*'/"^ to find, 

U'it,x) := e'(«-"+^*)/^ J e*((-(i+^C)— (i))*/e+C^) ni(T(^ + eC), ^ + ^C) G(C) dC 

:= e'(«-"+^*)/^ a(£,t,x). (5.1) 

Expanding in e and keeping just the leading term yields the principal term in the geometric optics 
approximation 



ije _ g»(i-^+rt)/s y e*("<-"(i)<*) nL(T,0 G(C) dC, v(0 :=-a5T(0. 

One has 

jje ^ ^Jii.+rt)/e ^^(^ „ ^(^)^) ^ J ^^.< J^^i^^Q G(C) dC . 

A complete Taylor expansion yields the corrected approximations which satisfy the equation with an 
error 0{e^) for all N. We write 0{£°°) for short. This yields infinitely accurate solutions, 

U^it,x) := e'(i'=+^*)/^a(i,a;,e), a{t,x,e) - ao{x - vt) + eai{t,x) + ■ ■ ■ . (5.2) 

If is a semi-simple eigenvalue of L(t, ^), and $o S Ker L(r, ^) \ {0} of dimension 1, then the leading 
amplitude gq in the case of (jl.ip (resp. (jl.Sp ) is of the form 

a(t,x)$o, (^resp. a(t, x) Ai$o, ■ • ■ , ^^rf^o^ 
with scalar valued amplitude a satisfying, 

{dt+w.da:)a = 0. 

This shows that a is constant on the rays which are lines with velocity equal to the group velocity 
v(0 := -d^riO. 

For g G \ the solutions do not have compact spatial support. This weakness is easily overcome. 
Choose X € C^(R'') with x = 1 on a neighborhood of the origin. For g G S{M.'^), define exact solutions 
by cutting off the integrand outside the domains of definition of t(^) and ni(T(^),^), 

u%t, x) :^ / e'(«-+-(«)*) n,.(r(0, - ^/e) x{V~e{^ ~ ^/s)) • (5.3) 



The analysis above applies with the only change being the initial values. In the preceding case these 
values were equal to the transform of ni(r(^),^) g{^ — ^/s) and in the present case they are infinitely 
close to that quantity, 

u%0,x) = J e"-« TIl{t{0,0 ai^-i/e) d^ + 0(£°°) • 

This yields infinitely accurate approximate solutions (2.5) which have support in the tube of rays with 
feet in the support of / e*^-« Ul{t{0^0 9(0 d^- 



5.2. Geometric optics with variable coefficients 

The Fourier transform method of the preceding sections is limited to problems with constant coefficients. 
In this section the WKB method which works for variable coefficients is introduced. It will also serve for 
the analysis of reflected waves. 

Let £ be the general operator in ()2.1|) . Fix (r, ^) G Char£ and seek asymptotic solutions 

- e*'^/^^e^aj (t,a;), with the phase ^(t, x, C) = tr + x^. (5.4) 

3=0 



Contents 4J 

May 9, 2011 5:49 WSPC/INSTRUCTION FILE confluentesvf 
More precisely we construct smooth functions aj{t,x) with suppa^ n ([0,T] x R'') compact so that if 

CJO 

a{t,x,£) ^ e-' aj{t,x), 

3=0 

in the sense of Taylor series at e = 0, and suppa n ([0,r] x M''x]0,e]) is compact, then 

:= e''^/^ a(<,a;,e) 

satisfies for all s,N, 

Wcu'W = o(£^) 

11"^^ llff=([0,T]xR<i) 

In this case we say that ()5.4p is a inGnitely accurate approximate solution. The next result recalls some 
facts about such solutions. 

Theorem 5.1. Suppose Problem (j2.1[) is hyperbolic, (r, ^) £ Char(>C) satisfies the smooth variety hy- 
pothesis, and is a semi-simple eigenvalue of C{t^^). 

(i) If the coefficients aj satisfy the recursion relation 

ao{t,x) e Ker £i(t, (5.5a) 

Vj>0, iCi{T,Oa,+i{t,x)+C{dt,d^)a,{t,x)=0, (5.5b) 

then (|5.4p is an infinitely accurate approximate solution of (|2.ip . 
(a) If Qj = n£(T, ^)(7j G C(J°(MJ^) are supported in a fixed compact K, then there is one and only 
one family of Uj satisfying (j5.5p together with the initial conditions, IIc{t, ^)aj{0, ■) = gj and the 
polarization Ilc{T,^)ao = ag. They have support in the tube of rays with feet in K and speed of 
propagation v(^) = — 9^r(^). 
(Hi) The principal term ao is a solution of the transport equation 

dtao + v(C) • 9,ao + Tlcir, OBix)IlciT, ^ao = 0. (5.6) 

Proof. For simplicity note He := IlciT,£,) when no ambiguity is to be feared. 

The equations ()5.5|) are obtained by injecting in ()2.ip . to find an expression ^ e''^/^ X]j>o ^"'^j 
In order that the Wj vanish it is necessary and sufficient that the equations (|5.5p are satisfied. 

Next examine the leading order terms to find the relations determining oq. Projecting the case j ~ 
of (|5.5p onto Ker£i yields, 

d 

He {dt+Y,-^i9i+B{x))ao^0. 
1=1 

This yields a first order system satisfied by oq — He oq, 

d 

dtan + ^ ^, dian + Uc B{x) lie ao = 0. (5.7) 

1=1 

The leading order part of this equation is a scalar transport operator. To see this differentiate 
>Ci(T(0,C)n£(r(0,C) = with respect to to find 

Ai + ^Id) nar(0,C) + A(r(0,C)^ (n£(r(e),0) = 0. 
Multiplying on the left by n£(r(^),^) eliminates the second term yielding, 

HcAiWc + ^^n^ = 0. 

Injecting this in (|5.7I) yields 
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In order to compute the coefficients recursively, muhiply (|5.5bp on the left by the partial inverse 
QciTjOi using the identity in (|2.7p . to obtain for j > 1, 

iI~Uc)aj = iQcCidt,d^)aj^i. (5.8) 

Projecting (j5.5bp on the kernel yields, 

Writing aj as 

a-j = Uca-j + (/ - Ilc)a-j, yields Uc C{dt,dx)Tlc a-j = -IlcC{dt,dx) {I - Uc) aj. 
This is again a transport equation, but with a righthand side, 

dtUcaj + V • dxUcaj +ncBUc a, = -Uc C{dt,dx) {I - Uc) a,. (5.9) 
(|5.8p and (|5.9p permit to calculate the coefhcients recursively, knowing the initial values. □ 



Next apply the above algorithm to the PML operator L. Fix (r, ^) e CharL and seek asymptotic 
solutions 

^ e''^'^^e^hj{t,x), with the phase S{t,x) = tr + x ■ (5.10) 

Corollary 5.1. Suppose Problem (jl.ip is strongly well posed, (t, ^) G CharL satisfies the smooth variety 
hypothesis, and is a semi-simple eigenvalue of L(t,^). 

(i) If the coefficients dj satisfy the recursion relation 

ao(t, x) e Ker Li(t, ^), (5.11a) 

Vj>0, {I~Uj^)dj{t,x)^iQiL{dt,dx)&j-i{t,x), (5.11b) 
dtUj^ aj + V • dxU^ hj + I3{x) % dj = -Hj L{dt, d^W - Uj) dj. (5.11c) 

then (|5.10p is an infinitely accurate approximate solution of ()1.5p . 

(ii) If gj{x) = Uj^ cjj £ C^{W^) are supported in a fixed compact K , then there is one and only one family 
of dj satisfying ()5.1ip together with the initial conditions, UL{T,^)dj{0,x) = gj and the polarization 
^LiT,i)dQ = flQ. They have support in the tube of rays with feet in K and speed of propagation 
v = -9«t(0. 

(Hi) The principal term do is a solution of the transport equation 

dtdo + v-dxdo + l3{xr^o = 0, wzth pix) = V^^%^^. (5.12) 

Tit] atl 



1=1 



Proof. We need only identify the constant term in (|5.6p . Use the form of the projector given in Propo- 
sition [27T1 to obtain 

U^^Bix)U^ ^Pix)U^. □ 

5.3. Amplifying layers 

The coefficient cri(a:i) > is introduced with the idea that waves will be damped in the layer. In this 
section, we show that sometimes the anticipated decay is not achieved, and waves may be amplified. This 
was observed in [5]. The authors analysed the phenomenon for a constant in the layer. They showed that 
in an infinite layer solutions can in certain cases grow infinitely large. We present a related analysis using 
WKB solutions which has three advantages, 

1. The analysis is valid for variable coefficients cri(xi) which corresponds to common practice. 
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2. The growth is seen immediately and not expressed in terms of large time asymptotics. 

3. The analysis in [5| was in part restricted to d = 2 and eigenvectors of multiplicity one. We remove 
these restrictions. 

It is because of 2 that wc choose not to follow the authors of [5] in calling this phenomenon instability. 



Theorem 5.2. Suppose (r, ^) e Char(L) satisfies the smooth variety hypothesis and j3{x) is as in (|5.12p . 
Suppose in addition there is an interval on a ray 

T := {{0,x) +t{l,~d^T{0), 0<t<to}, sothat, [ " Pix- td^ri^))) dt < 0. 

Jo 

Then the corresponding WKB solution grows in the layer. 

Proof. The solution of the transport equation (|5.12p evaluated on F is 



ao(io,a;) = exp j f3{x- sd^r {(,)))) ds j ao{0, x + tod^ri^)). 

The exponential is strictly greater than 1, so 

\ao{to,x + tod^T{0)\ > \ao{0,x)\. 



□ 



Example 5.1 (No amplification for Maxwell/D'Alembert). If the dispersion relation is = 
and CT > then there is no amplification since, 



13 



d ^2 



i=i 



ll^li^ 



> 0. 



Example 5.2 (Amplification is common). For the dispersion relation = q{£^) where g is a positive 
definite quadratic form so that the ^ axes are not major and minor axes of the ellipse q = 1, there are 
always t > 0,^ so that xi layers with cri > arc amplifying (jU). There are two lines on {r ~ qiO^^^} 
where dq/d£,i = 0. The half cone on which dq/dS^i < corresponds to rays on which xi is increasing so 
they enter a layer xi > 0. The half cone {dq/d£,i < 0} is divided into two sectors by the plane = 0. 
The sector on which > (resp. < 0) corresponds to growing (resp. decaying) solutions (see Figure 
[1] on the left). This example shows that amplification is very common. Consequently for the dispersion 



outgoing, ^ < 





Fig. 1. Amplified outgoing wave numbers in bold 



relation ~ q(^) it is wise to align coordinates along the major and minor axes of the ellipse to avoid 
amplification. However, if (r^ — qi){T'^ — 92) divides the characteristic polynomial and the axes of qi and 
q2 are distinct from each other then no linear change of coordinates can avoid amplification in the layer. 

A second example from |10| is the linearized compressible Euler equation with nonzero background 
velocity (c, 0), c > for which amplified wave numbers at a right hand boundary are indicated in bold in 
the right hand figure [TJ 
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Summary. There is no amplification when the characteristic polynomial is a product of factors r and 
— q where g is a positive definite quadratic form with axes of inertia parallel to the coordinate axes. This 
includes the cases of Maxwell's equations in vacuum, for which the method was developed by Berenger, 
the linearized Euler equations about the stationary state, and the linear isotrope elasticity equations. For 
these the quadratic forms q are multiples of 

Example 5.3 (Methods related to Berenger, continued). For the model developed in Section [321 
for the Maxwell equations, one can compute 

/? = 2^ ^ a,ix,)^ > 0. 

Thus, this model has exactly the same good properties as Berenger's, and is strongly well posed. For 
Maxwell equations, it is therefore an attractive alternative. The advantage is two fold. The system with 
the auxiliary variables is very compact. And it is strongly well posed, even for discontinuous a. 

On the surface this result sounds almost too good to be true. However the Berenger system in the case 
of Maxwell's equations has almost exactly the same structure. The energy method proof when cr" S L°° 
shows that there is a large vector V consisting of the components of U together with differential operators 
Pa{D) applied to U and a strongly well posed equation for V. This means that if one were to introduce 
the additional variables in V one obtains a system with some of the desirable properties of SPML (strong 
PML). However, the SPML reduction is much more compact, and, has a good energy estimate even 
when a is discontinuous. The extension of this strategy to other equations is not straightforward. For 
elastodynamic models, see |28) . 

6. Harmoniously matched layers 

This section introduces a new absorbing layer method. It is based on the following strategy. Start with 
an operator L = Li{d) on the left and consider a smart layer on the right 

R{t,x,d)=Li{d) + C(t,x), C ^ a{xi)[T:+{Ai) + vTT^{Ai)^, supper C [0, oo[, (6.1) 

generalizing (jl.2l) . This method is embedded in a family of absorbing layers parameterized by /i > 0, 

i?^ Li + mC- (6.2) 

The method is nonreflective when /i = and is both reflective and dissipative for ^ > 0. When a 
is discontinuous the leading order reflection coefficient for wave packets of amplitude 1 oscillating as 
^i{Tt+xC)/e jg q£ ^j-^g form e^r(T, ^). The leading order reflections can be removed by an extrapolation 
method using two values of ^. This simultaneously removes the leading reflections at all angles of incidence. 
We call the resulting method the harmoniously matched layer. 

6.1. Reflection is linear in /j, by scaling 

In this subsection the linearity in /i of leading order reflections by the layer with from (j6.2[) is demon- 
strated by a scaling argument when a{xi) = l^j^^o- In the next subsection the reflection is computed 
exactly for Maxwell's equations yielding additional information. 
If (Li + C)U = then 

U{t,x) --Uifit^fix), satisfies i?^C/ = 0. 

Suppose that U has an incoming wave of wavelength e and reflected waves Ui with amplitudes piS. 
Then U_ has an incoming wave with wavelength e :~ £/ p- The reflected waves have amplitudes 

£ 

PiE = pep - = PiP £■ 
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Denote by the reflection coeflicient of R^. The leading amphtude of the reflected i wave is then p^e. 
The preceding identity shows that p^ = pgp showing that the reflection coefficients are hnear in //.El 

6.2. Reflection for Maxwell with smart layers 

In this section £ may denote one of two distinct operators. One option is the Maxwefl operator Li from 
p.2|) for the C'^ valued field E + iB. The lower order term is B := p-C from the smart layer ()6.1|) . 
Alternatively C may denote the Berenger operator operator L with lower order term B = pC with 

^a(xi)/jv ... 0\ 
C = I ; '■. : ' suppo- c [0,oo[. 
... 0/ 

In both cases the absorption term is linear in p. We compute the dependence of the reflection coefficient 
on p. 

Lemma 14.61 shows that the Cauchy problem is equivalent to homogeneous problems in each half-space 
with a transmission condition on F := {xi = 0}, 

[AiU]^ = 0. (6.3) 

In order to cover both cases the operator, coefficients, and unknown are indicated with round letters. 

We study the reflection of high frequency waves in xi < which approach the boundary xi — 0. The 
input is an incident wave with phase S^{t,x) :~ rt + £^'x, where r 7^ and r(^) = ±|C|. The phase is 
chosen so that the group velocity v = — ^/r satisfies vi > 0. Denote x' := {x2,X3), — (^2,^3). Theorem 
15.11 applies to the incident wave with i3 = 0. Ina;i<0, the incident wave is 

CXD 

W ■.= e''^'^*''''>/'a'{t,x,e), a'(i, a;, e) - ^ e^' aj(t, x), /:i{dt,d.^)W ^ 0{e°°) . (6.4) 

3=0 

Suppose that the amplitudes arc supported in a tube, T, of rays with compact temporal crosscctions 
rn {t = 0} CC {xi < 0}. 

We construct a transmitted wave with the same phase, and a reflected wave with phase S'^itjx) := 
rt + ^'^x, with ^'^ := {—£,!,(,')■ Wc first show that there are uniquely determined reflected and transmitted 
waves. Then we compute exactly the leading terms in their asymptotic expansions. 

The reflected wave is also supported in xi < 0. The group velocity for the reflected wave is equal 
to v" :~ (— Wi, v'), and in xi < 0, 



= e*^''(*'")/^a«(<,x,e), a''{t,x,e) ^ af{t,x), Ci{dt,d^)V' = 0{e^) . (6.5) 

3=0 

The transmitted wave is supported in xi > 0, 

CXD 

:=e'^'(*^"'/"a^(<,x,e), a^(<,x,e) - ^ £^'aj(t,x), C{dt,d.,)W' = 0{e^) . (6.6) 



3=0 



Theorem 6.1. 



(i) Given the incoming amplitudes a'j there are uniquely determined amplitudes aj and so that 
for any choice of the a^'^''^{t,x,e) ~ ci''^'^it,x), theU^, and are infinitely accurate 

^This argument can be made rigorous under the following conditions. The incoming wave is a wave packet with oscillatory 
part e'('^*+^^'/'^ with (t, ^) G CharL. Denote (r, the part determining the oscillations in xi = 0. Consider the roots of 
each of the equation, det ii (t, , 5') = 0. The nonreal roots are called elliptic. They lead to waves which have the structure 
of a boundary layer of thickness ~ e. The real roots are called hyperbolic. The favorable situation is when all the hyperbolic 
roots are at smooth points of the characteristic variety and the group velocities are transverse to the boundary. In that case 
one can construct infinitely accurate asymptotic solutions of the transmission problem consisiting of incoming, reflected, 
and transmitted wave packets corresponding to the hyperbolic roots, and, a finite number of boundary layers corresponding 
to elliptic roots. As this is a long story, we content ourselves with the Maxwell computation of the next subsection. 
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solutions of the differential equations and the transmission condition is also satisfied to infinite 
order, 

\f{t, x') eRxR^, Ai{W +V'){t,0_,x') = Ai{W'){t,0+,x') + O{e°"). (6.7) 

(ii) In the case of Berenger's PML, the coefficients vanish identically for j > 0. 

(iii) For the smart layer i6.1\ \6.2^) with a = X^^^q, the coefficient vanishes identically. The reflec- 
tion coefficient of the layer is equal to 

, ^ £? — T'^ u(l + I') ?j? — 1 



8T£_f 8r vf 

That is, z/a5(t,0_,.T') = a{t, x')<^{t, G Ker L(r,C), then 

a''{t,0_,x') ^ ei?(r,e)a(t,a;')$(r,C«) + 0{e^) . 

Furthermore, the amplitudes a^ " are such that on the interface T, we have for all G N^, 

i = and j < 1, or i > 1 and j > 00 

dia[~dial G fi{C,+,.l[^J.](g>C^), diaf G m(Q+j-i[m]®C^)- (6.8) 

(iv) The smart layer with a{xi) satisfying ct(0) = • • • = crC^^-'^) (0) = 0, a^''^ (0) ^ is nonreflecting at 
order k for any angle of incidence, i.e. if aQ{t,0_, x') = Q;(t, a;')$(r, ^), there exists Rkir,^) such 
that 

a«(i,0_,x') = £'^(TW(0)i?,(r,e)a(<,a;')$(^,r) + 0{e''+') . 

Furthermore the amplitudes a^' " are linear functions of ^ on the interface T. That is denoting 
cf (^) = af |p and cf{fi) = aj ^-p — al |p, we have for all i > k in N, 

cf'^ifi) G /i(Q_i[M]®C3). 

Remark 6.1. 1. There exist choices of a^'^''^ so that W^, V^, and is an exact solution. Since the 
transmission problem is well posed, there is a uniquely determined corrector c/ smooth and infinitely 
small on both sides so that adding c*^ yields an exact solution. Adding to the left corresponds to 
adding the infinitely small term c^ e*"^ to with a similar remark on the right. 

2. Part (iv) of the theorem with fc = generalizes part (iii) to discontinuous and variable cr(xi). 

3. The basis elements, for af and ^' for q'q are homogeneous of degree 2 in r, Doubling t, ^ and also 
£ leaves the incoming and reflected waves unchanged. Therefore ei?(r, <^) must be equal to 2£i?(2r, 2^). 
This explains why R is homogeneous of degree -1. 

4. The reflection coefficient vanishes when ^' = 0. Since it is an even function of V^i? = too. 
Proof. The incoming solution is given. 

(i) Seek the leading amplitudes ag and Oq. We will show that Oq = so it is actually af that is the 
leading amplitude of the reflected wave. A jump discontinuity in a lower order coefficient does not lead 
to reflection at leading order. Denote 

Ct ■■= dt + A2d2 + Asds ; Ci Ct + Aidi ; C Ci + fiC . 

T '.= dt +V2d2 +V3d3 is the tangential transport operator. By Theorem l5.ll the amplitudes are polarized, 
i.e. Oq ^ = n^ag ^, and Og (resp. ag ) satisfies a forward transport equation in xi > (resp. backward 
in xi < 0) with zero initial values in time, 

(-fiai+r)af = 0, xiGM-, (6.9) 

{vidi +T + fiUcCUc)a^ ^0, xi G M+. 



'^Cj [p] denotes the space of polynomials of degree less than or equal to j with complex coefficients. Cj [fi] (X) is the 
corresponding space of polynomials with coefficients in C'^. 
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Therefore, to determine Og and Og everywhere, it suffices to know Cg (i, 0+, x') and Cg (i, 0_ , x'). These 
values are determined from the transmission condition (16.31): 



Ai{a'^{t,{)_,x') + a^{t,{)_,x')) ^ Aial{t,{)+,x'). (6.10) 

The matrix Ai is singular. It is easy to see that (Ker£(T, ^) © Ker£(r, ^^)) n Ker„4i = 0. Therefore, 
^iKer £(t, ^) and ^iKer £(r, are complementary subspaces and generate Range A\ . This proves that 

a^{t,Q_,x') = 0, al{t,Q+,x') - a^(f,0_,x') = 0. (6.11) 

By the transport equation, we conclude that 

Og = 0, for xi < 0. (6.12) 

The reflected zeroth order term vanishes identically when x G Mil. We also deduce from the transport 
equation (|6.9p that 

vi{dial ^ dia'o) + /iHiCnLaJ = on T. (6.13) 

Next determine inductively the correctors. For simplicity, throughout the proof we note n£ := Hc{t, ^) 
and := ^l{t, C"). Write the recursion relation (|5.5I) for j >1 for the incident, reflected and transmitted 
waves. Split the amplitudes as 

aj-'^{t,x) = n£a^'^(i,a;) + [I - He) a'j'^ (t, x), 
af{t,x) = Il"caf{t,x) + (I-U^) af{t,x). 

(I — He) ajit, x) and (/ — H^) a" (i, x) are determined directly by (|5.8p . To determine the projection on 
the kernel, split the transmission condition (|6.3p and insert (|5.8p on the interface to get, 

Ai{llca'j{t,0,x') ~Uca^it,0^,x') + n^afit,0_,x')) = 

- Ai{{i - n^) a'j{t, 0, x') - (/ - n^) aj{t, 0+, .t') + (/ - n2) af{t, 0_ , x')). 

As for the terms of order 0, this determines Hcajit, 0+, x') and Il^a^{t, 0_,x'). By (|5.9[) . the projections 
are solution of a transport equation, therefore uniquely determined by initial data and the values on the 
boundary . Borel's theorem allows one to construct 

a'(t,x,e), a^(t,x,s), and, a"(i,a;,e), 

so that the transmission condition is exactly satisfied. With this choice the approximate solution satisfies 
the transmission problem with infinitely small residual. 

(a) Theorem l4.4l implies that the exact solution in xi < is equal to + 0{e°°). The error of the approxi- 
mation is 0(e°°) so the exact solutionis equal to Z^^ + V-+0(£°°). Therefore V = {W + V')-U' = 0(e°°) 
which is the desired conclusion. 

(iii) For the smart layer (|6.1[ 16. 2|) with a = \xx>Oi compute the first order term by (jS.lip with j ~ 1. 
First deduce from (|5.8p that 

[I -T\L)a{{t,x) = iQLLi{dt,d^)a'oit,x), xi e R, 

(/-n£)af(i,x) = 0, xiGK-, (6.15) 

(I -UL)aUt,x) ^ iQLiLi{dt,d,)+fiC)a^it,x), xi G M+. 

Replace in (j6.15p the xi derivatives using ()6.9p . 

{I-nL)a{{t,x)=iQL{LT + Aidi)a'^it,x)=iQL{LT~^Air)a'o{t,x), 

{I - UL)alit, x) = iQL{LT + Aidi + M C) a^{t, x) = IQl {Lt - ^^AiT + ^i (-^^iHl CHl + C)) al{t, x) , 
to obtain, with Ci := C - — AiIIl CIIl, 

Vl 

(/-ni)a[(t,0+,x') - (/-nL)ai(t,0,a;') = i^iQLCia'S,Q,x'). (6.16) 
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Using (/ — Ili^)af = in the transmission condition yields, 

Ai(n«af + ULa[-IlLal) = ifxAiQLCia'^ . 



(6.17) 



The eigenvalues of Ai are and ±1, with associated ortlionornial set of eigenvectors $o = and 
$± = (0, 1, ±i)/-\/2. The projection operators on the positive and negative eigenspaces are 7r±(Ai) = 
^±^±1 and C = ^+^*_^ + z/4>_$l. The kernel of L{t,^) is one-dimensional, it is spanned by 



(6.18) 



and the projection on KerL(T, ^) is 11^ = Compute 

<f><f>* <f><f>* 



^2-$($*4'+)($!^$)** + ;/$($*$+)($;$)$* 



($*$)2 



Define 



Since Oq is polarized. 



so, UlCUl = . 



(6.19) 



Cia'o = {C- ^^AiUl) a^, = (C - ^Ai) := Ci a^. 



To compute the righthand side of (|6.17p . use 
^0 \ 



I 

2 2 

1/ - 1 1/ + 1 



,7 = 



(r-ei)^ + ^(r + 6)^ 
4t2 



-((l+t,l)'+Kl-^'l)'): 



to find 



Ci = 



v + 1 










1 

I'l + l 

2vi 



\ 

.vf + 1 

I 

2wi 

1 



Write = oi^'^, and compute 



Ci $ = (1 + i.; 



to find a new version of (j6.17p . 



e2 _ 2 

Ai(n£af + n^ai - nLa[) = + z.) ^^Vt^ "o ^iOl * ■ 



(6.20) 



Next compute Ql ■ First compute a basis of eigenvectors for L{t, ^). $2 is such that L(t, ^)<I'2 = ^$2, 
$3 is such that L(t, ^)$3 = 2t$3. Choose 



Note that 



^2^L $3 = 3'(-T,0- 
* = re-ei(Tei+ieAei), 



which gives 



and 
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and 

$ = ^ - -^(rei -i^A ei), <i>3 = ^ - ^{rei + i£, A ei), 
?i SI 

(rei+z^Aei) = -(^-$3), 

T 

T T r 

Since Ql is the left inverse of i, we have Ql^ ~ and Ql'I's = ^'^'sj which gives 

2 _ pi p2 
Ql^ = + ^,$3. 

Write the coefficients on T as 

TlLa[^ = a[^$, n«af = af , 
and inject into the transmission condition to obtain 

Since the kernel of Ai is ei, Ai in the preceding identity may be replaced by the projection on (e2, 63). 
The projection of $ is = ^' — i-^^ A ei, and note that $3 and $(r, have the same projection, which 
is (p3 = + i-^S, A ei. Write 



a1^3 + [a{ - al)^ = + v) (^2^?' + T^^Vs) 



- - 1 



The solutions are parameterized by Uq, 

^2_^2 /^2_^2-)2 

Now use the results in (i) and prove (j6.8p by induction on i. Equation (|6.12p asserts that ag = for 
xi < 0. Equations (|6.9p . and (|6.10p . imply that at xi = 0, = IlLaQ^, and 

vi{diaQ - dia'„) + T{al - «„) + /i7(ao - &„) + ^700 = 0, xi> 0. 

Differentiation in xx several times yields d{aQ — d{aQ E fiCj^ilfi] ^ C^, giving the results for i = and 
J > 1. 

Assuming the inductive hypothesis is true for i we prove it for i + 1. Write (jS.lip in the form 

(/ - nf ) af+i(<, x) = iQl{LT{dt,d,,) + Aidi) aiit, x), 
{I-Hl) af+i(t,a;) = iQULridt^d,,) + Aidi) al{t, x), 
{I-Hl) ar+i(t,x) = iQL{LT{dud^,) + ^i^i + ^iC) aj{t, x). 

By induction, dl{I — Hf) af{t, x) E /iCi+j_i[(T] (8) C'^ on the interface. Write for xi > 0, 

(/ - Hl) (af+i - a^+i)(t, x) = iQL (Lt + A^di + ^iC) (af - a]){t, x) + ttiQLC {t, x). (6.21) 

The inductive hypothesis, shows that on F, d{{LT + Aidi+ iiC) {aj — a\) E /iCi+j[/i]®C"^ and d{{^iC al) G 
/i Co [cr] ® C^. The result follows for (/ - Hl) ( a,J_^_l — al^i), and for (/ — n£) Q^f+i hi the same way. The 
transmission condition extends the assertion to the other parts 11^ {cij_^i — al_^i) and n£ af_|_j^. 
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(iv) Here a vanishes to order fc at xi = 0. (|6.1ip and (|6.12p are still valid, and the transport equations 
6.9p implies on the interface F that 

dinLa^ - diULa'o = 0, = 0, • • • , fc, 



From (|63T|) for i = 0, (|05l) and (|02|) . derive 



af = Xliaf everywhere, 

5^(/-ni)a^-c)J(/-ni)a{ =0, j = 0,---,fc-l, on F, (6.23) 
9f (/ - Hi) a[ - (/ - Hl) a{ = ina(''HO)QLCiIlLa'o on F. 

Using the transmission conditions to obtain on the interface F, 

n£af = 0, ULa[ = Uia^. 
Insert into the transport equations (|5.9p to find = in R_ , and 
d(llLal - dlllLai = 0, j = 0, • • • , fc - 1, on F, 

afn^a^ - d^ULa{ = - J-niAi(9f(/ - n^) al - d^I - n^) a{) = -iiifT«(0)ni^iQLCinLa^ on F. 

(6.24) 

Recover 

d{al - d{a{ = 0, J = 0, • • • , fc - 1, on F, 

- d'la{ - i/xaW(0)(/ - ^nL^i)niAiQiCinia^ on F. 

Now proceed iteratively, to sec that, for z < A; + 1, 

af = n^af = everywhere, 

d(af-dia{ = 0, j = 0, on F, 

d1-^+\l-UL)aJ~d1-'+\l-UL)ai=iQLAiidt'+^aJ^^-d'^-'+^aU) on F , (6-26) 
a^'+iHiaf - dt'^^ULai = -l-ULAi{dt'+\l - Ul) aj - a^*+i(/ - Ul) a^) . 

Denote by Si the value of cJf ^'^^af — cjf ^'^^a^ on F. Equation (|6.26p yields the recursion relation 

s, = i{i - — nL^i)QLAiSj_i , Si = i^lCJ^^\Q){I - — niAi)QiCinia^ , 

Vl Vl 

which can be solved as 

Sk+i - i^i(J'^''\f^)^mLal with M (i(/ - ^I{LAi)QLAifQLCi. 
The first nonzero refiected term is therefore a^^i = Wf^a^^^, and using the transmission condition yields 

v4i(nLa^+i + nL4+i-nLa^+i + sfc+i) = o. 

The incoming amplitude on F is Aq = Q;(i,x')$, the leading reflection is a^^j^ = a^^j^^" and the leading 
transmission is a^^j^ = q;^_,_j^$^. Using again the notation $' to denote the projection of a vector (f> on 
Vec{e2,ez), this linear system is solved as 



a 



'^+1 $'A($^)'' "''^^ $'A('I>«)' ■ 



■ (fc).nN *'A(M$)' ^ , . ffc.,^, ($«)' A (A/r$)' 

K+i A ($^)' ^ <&' A ($"■)' 

The proof of the linearity follows the same path as in {Hi). □ 
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6.3. Harmoniously matched layers 

Based on Theorem 16.11 we construct an extrapolation method for symmetric hyperbolic operators with 
smart layers which eliminates the leading order reflection. The resulting method has desirable stability 
properties and is nearly as good as Bercnger's algorithm for the Maxwell equations where his method is 
at its best. We think that the new method provides a good alternative in situations where Bcrenger's 
method is not so effective. 

Consider the computational domain xi < bi. The domain of interest is the interval xi < ai < bi. 
The absorbing layer is located in ai < xi < 6i. The differential operator in the computational domain is 
symmetric hyperbolic L with smart layer 

LU + (Ti{xi){tt+{Ai) + i^TT^{Ai)) U ^ , cti > , supp cti C {xi > ai} . 

At the outer boundary xi = bi of the absorbing layer impose the simplest weakly reflecting boundary 
condition 

TT-{Ai)U = when xi = bi . 

This is a well posed problem provided that Ai has constant rank on xi = bi. When L = Li{d) has 
constant coefficients it generates a contraction group in L'^{{xi < bi}). 

The hamoniously matched layer algorithms compute a smart layer with coefficient <ti and also 
with coefficient 2ai. In view of Theorem 16. 1[ subtracting the second from twice the first, 2U{ai) — 
U{2ai), yields a field with one more vanishing term in the refiected wave at the interface xi = ai. This 
extrapolation removes the leading reflection. 

The harmonious matched layer algorithms in a rectangular domain TZ perform the same extrapolation 
with absorptions in all directions. With 

d 

LU + ^ (Tj(a;j)(7r+(Aj) + vTr^{Aj))U — 0, <Tj > 0, supper^ C {\xj\ > aj} . 

with 

7rzp(j4j)C/ ~ when Xj — ±bj . 

This initial boundary value problem on a rectangle has weak solutions. El When L = Li{d), the L^TZ) 
norm is nonincreasing in time. The extrapolation is 2 U{ai, • • • , cjd) — U{2(Ji, • • • , 2ad)- 

Open Problem. For discontinuous Cj, the uniqueness of solutions to the initial boundary value problem 
on the rectangular computational domain is not known because of the discontinuity of the boundary space 
ker Aj at the corner. Solutions are typically discontinuous. Uniqueness of strong solutions and existence of 
weak solutions is proved by the energy method. We do not know how to prove uniqueness of solutions with 
regularity not exceeding that of solutions known to exist. Similar problems plague virtually all methods 
on rectangular domains with absorbing boundary conditions imposed on the computational domain with 
corners. The present problem is one of the simplest of its kind. The fact that algorithms designed to 
compute solutions encounter no difficulties is reason for optimism. 

6.4. Numerical experiments 

Simulations arc performed for the 2-D transverse electric Maxwell system in the (a;, y) coordinates, 

dtE^ - dyH, = , 

dtEy + d,H, - , (6.27) 

dtH, + d^Ey - dyE^ = , 

"^This can be proved by penalisation. Denote by Q the rectangular computational domain. Add Al^d.^^ to L and solve on 
The limit as A -> oo provides a solution in L°°([0, T] : L'^ifi)) [i]. 
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in a rectangle, with boundary conditions n A E = on the west, north and south boundaries. The layer 
will be imposed on the east boundary. Maxwell Berenger is given by 



dtE, - dyH, = , 

dtEy + d^H, + a{x)Ey = , 

dtH,, + d^Ey + a{x)H,^ = , (6.28) 
dtHzy — dyEx ~ , 

For the computation, these equations are used in the whole rectangle (see the discussion in the introduc- 
tion), with cr = outside the layer. The boundary conditions are 

Ey = Hz and E'^ = on the east , n A -E = on the other boundaries. (6.29) 

Since n+(Ai) = (0, 1, 1) and n_(^i) = (0, 1, -1), the smart layers are: 

dtE., - dyHz = , 

dtEy + d,Hz + ^{Ey + Hz + v{Ey - H ,)) = , (6.30) 
dtH, - dyE, + d,Ey + ^{Ey + Hz- v{Ey - Hz)) = . 

The boundary conditions (|6.29p are imposed. 
The Yee scheme for Maxwell is 



At Ay 



- = 



At Ax 



A< Ax Ay 

The Yee scheme for Maxwell Berenger using the notations cTi = (7{xi) and cr.i^i = cr(xj_,.i) is. 



= 0. 



(£^x)" 1 - - (£^x)T' ■ (Hz)":? 1 - {Hz)"":? . 1 

At Ay ~ ' 

At + Ax ^ 2 



0, 



= + f^^+i n = 0: 



At Ax 



At Ay ' 
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The Yee scheme for the smart layers is 

M Ay 
At Ax 

^ (i+.k(^.)".+i+(^)";:^^ (i-.)^.+^(^)::;^,^i+^.-i(^):r7,^. 

2 2 2 2 
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(6.33a) 



, (6.33b) 



At Ax Ay 



= 0. 

(6.33c) 



The schemes are implemented using time windows to save memory. 

The harmoniously matched layers can be implemented in several ways that we compare. The function 
cr{x) is as above. 

HML Version 1. Global extrapolation. Compute the solution of (16.33^ with an absorption of a, 
iE\H'^) and 2a, [E'^.H'^) over the whole time window. Then E^^y = 2 * El y - El y and 
H-. = 2* - H'^. 

z z 

HML Version 2. Local extrapolation. Compute at each time step the solution of (|6.33p with an 
absorption of cr, {E^,H^) and 2a, {E^,H^) over the whole time interval. Then E^.y = 2 * 
E], y — El y and Hz ^ 2 * Hi — H^ . Save computation by taking advantage of the fact that the 
computation of E^ does not involve the absorption parameter. At each time step, 

(1) Ex is computed by (|6.33ap . 

(2) two values of Ey are computed by (|6.33bp : Ey with an absorption parameter equal to a, Ey 
with an absorption parameter equal to 2a. 

(3) two values of Hz are computed by ()6.33cp : Hi with an absorption parameter equal to cr, 
H^ with an absorption parameter equal to 2a. 

Then £:„ = 2 * E}, - Ef, and Hz = 2* Hi- H^. 

HML Version 3. Split field local extrapolation. At each time step, 

(1) -Ea; is computed by (|6.33a|) . 

(2) two values of Ey are computed by (|6.33bp : Ey with an absorption parameter equal to a, Ey 
with an absorption parameter equal to 2a. Then Ey ~ 2 * Ey — Ey. 

(3) two values of Hz are computed by (|6.33cp : Hi with an absorption parameter equal to a, 
H^ with an absorption parameter equal to 2a. Then Hz — 2 * Hi — H^. 

We perform a series of experiments to illustrate the transmission properties of the layers. The coef- 
ficient V is meant to achieve backward absorption and is taken equal to zero. The domain of interest is 
(0, 6) X (0, 10), the coefficient a{x) is supported in 6 < x < 10. The time of computation is 4, the initial 
electric field is zero. The initial transverse magnetic field, 

H^ = cos (tt ) cos(fc7rt; ) x\x-x,\<r 

is compactly supported in the ball B{Xc, r), with x^ = (5, 5) and r = 0.8. 

The time of computation is fixed such that there is no reflection on the exterior walls. The initial 
mesh is taken to be Ax ~ Ay ~ 0.1, At — 0.0702, and then divided by 2 twice. 
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In the first set of experiments, the absorption coefficient is constant in the layer, equal to 2. The initial 
magnetic field hits the layer at incidence 0°{ v = (1, 0)) or 45° { v = {1, 1)). 

In Table [T] we compare the performances on a high frequency wave (fc = 10), while in Table [2] we 
consider a low frequency wave (fc = 1). 

In Table [T] we compare the performances on a high frequency wave (k = 10), while in Table [2] we 
consider a low frequency wave (fc = 1). 



Table 1. Comparison of the L°° errors for high frequency, discontinuous 
absorption. 





normal incidence 


45° incidence 


refinement 





1 


2 





1 


2 


Berenger 


9.4e-02 


3.9e-02 


7.9e-03 


1.3e-01 


2.9e-02 


5.6e-03 


Smart 


5.2e-02 


1.3e-02 


5.1e-04 


6.2e-02 


l.le-02 


5.3e-03 


HMLVl 


3.4e-02 


3.1e-03 


2.1e-05 


4.5e-02 


1.2e-03 


5.5e-04 


HMLV2 


2.5e-02 


6.0e-03 


1.2C-03 


7.4e-02 


l.le-02 


1.7e-03 


HMLV3 


2.1e-02 


4.2e-03 


5.1e-04 


4.5e-02 


5.3e-03 


5.7e-04 



Table 2. Comparison of the L°° errors for low frequency, discontinuous 
absorption. 





normal incidence 


45° incidence 


refinement 





1 


2 





1 


2 


Berenger 


1.5e-02 


7.1e-03 


3.5e-03 


1.3e-02 


6.1e-03 


3.0e-03 


Smart 


2.0e-02 


2.0e-02 


2.01e-02 


4.3e-02 


4.2e-02 


4.2e-02 


HMLVl 


1.7e-02 


1.60e-02 


1.6e-02 


3.4e-02 


3.3e-02 


3.2e-02 


HMLV2 


1.8e-02 


l.le-02 


6.7e-03 


3.1e-02 


1.9e-02 


l.le-02 


HMLV3 


4.3e-03 


2.6e-03 


1.4e-03 


8.2e-03 


4.8e-03 


2.6e-03 



In Tables [3] and SI we perform the same set of experiments, but the absorption coefficient is now a 
third degree polynomial in the layer, equal to {x — 6)'^/8. 



Table 3. Comparison of the L°° errors for high frequency, continuous absorption. 





normal incidence 


45° incidence 


refinement 





1 


2 





1 


2 


Berenger 


3.8e-05 


1.9e-07 


2.1e-09 


2.0e-04 


9.1e-07 


1.6e-09 


Smart 


2.7e-05 


2.2e-07 


1.7e-07 


1.7e-04 


9.0e-07 


3.1e-08 


HMLVl 


5.5e-07 


6.0e-08 


5.6e-08 


5.6e-06 


1.2e-08 


4.7e-09 


HMLV2 


6.8e-07 


6.5e-08 


3.1e-08 


2.6e-06 


8.1e-09 


2.8e-09 


HMLV3 


5.8e-08 


2.4e-09 


l.le-09 


1.5e-06 


9.5e-10 


9.0e-ll 



Table 4. Comparison of the L°° errors for low frequency, continuous ab- 
sorption. 





normal incidence 


45° incidence 


refinement 





1 


2 





1 


2 


Berenger 


6.2e-07 


3.2e-08 


7.8e-010 


5.2e-07 


2.9e-08 


6.5e-010 


Smart 


5.3e-04 


5.3e-04 


5.2e-04 


3.9e-04 


3.8e-04 


3.7e-04 


HMLVl 


1.6e-04 


1.6e-04 


1.5e-04 


8.6e-05 


8.3e-05 


8.2e-05 


HMLV2 


4.1e-04 


2.0e-04 


9.6e-05 


2.0e-04 


9.8e-05 


4.8e-05 


HMLV3 


l.le-05 


5.4e-06 


2.7e-06 


5.9e-06 


2.9e-06 


1.4e-06 
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The Berenger layer performs well on every frequency and every angle of incidence. Among the 3 
versions for the HML, the third version is the best, which should be analyzed thoroughly. 

Next compare the method on a gaussian initial value, supported in (0,6) x (0,10). Table [5] uses a 
constant absorption in the layer, while Table [6] uses the same smooth absorption as before. 



Table 5. Comparison of the L°° errors for a gaussian 
initial magnetic field, constant absorption. 



refinement 





1 


2 


Berenger 


1.5e-02 


6.7e-03 


3.3 e-03 


Smart 


3.4 c-02 


3.4e-02 


3.3C-02 


HML VI 


3.0e-02 


2.9e-02 


2.8e-02 


HML V2 


3.6e-02 


2.5e-02 


1.6C-02 


HML V3 


l.Oe-02 


6.6e-03 


3.9e-03 



Table 6. Comparison of the L°° errors for a gaus- 
sian initial magnetic field, continuous absorption. 



refinement 





1 


2 


Berenger 


7.5e-07 


2.0e-08 


8.3e-10 


Smart 


4.3e-04 


4.2e-04 


4.1e-04 


HMLVl 


1.3e-04 


1.2e-04 


1.2e-04 


HMLV2 


3.0e-04 


1.5e-04 


7.3e-05 


HMLV3 


8.8e-06 


4.3e-06 


2.1e-06 



Finally, take unstructured random initial value, supported in the ball centered at (5, 5) and of radius 
1. In Table [71 the absorption coefficient is constant in the layer, equal to 3. 



Table 7. Comparison of the L°° errors for a random 
initial magnetic field, constant absorption. 



refinement 





1 


2 


Berenger 


5.7e-02 


4.9e-02 


4.4e-02 


Smart 


6.7 e-02 


6.3e-02 


5.4e-02 


HML VI 


5.1 e-02 


4.5e-02 


4.0e-02 


HML V2 


6.4e-02 


3.0e-02 


1.9e-02 


HML V3 


3.2e-02 


1.5e-02 


6.7e-03 



In Table [H the absorption coefficient is a function of x in the layer, equal to {x — 6)'^/8. 



Table 8. Comparison of the L°° errors for a random 
initial magnetic field, continuous absorption. 



refinement 





1 


2 


Berenger 


l.le-04 


5.0e-05 


4.4e-06 


Smart 


7.2e-04 


6.9e-04 


6.4e-04 


HMLVl 


2.1e-04 


2.2e-04 


2.0e-04 


HMLV2 


5.0e-04 


2.7e-04 


1.2e-04 


HMLV3 


1.5e-05 


7.9e-06 


3.7e-06 



Summary. When comparing the reflection properties, the harmoniously matched layer, version 3, is 
competitive with the Berenger layer. For very regular data, the Berenger layers outperform everything. 
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The performance of the HMLV3 gives hope the method with its stronger well posedness, more robust 
absorption, and small reflection at all angles will be a good method where Berenger has proven less good. 
For example, for non constant coefficients and nonlinear problems. We have taken pains to make the 
comparison where Berenger is at its best. In 2D with a layer in a single direction the HML has an extra 
cost. Since there are 5 quantities to compute at each time step instead of 4 for Berenger. This is no longer 
the case in three dimensions, since both strategics have to split 6 unknowns. 

Open problems. 1. Our analysis does not explain the much better behavior with continuous absorption, 
nor the advantages of IIMLV3. 2. A comparison with other methods where only supplementary ordinary 
differential equations are added should be made. 
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